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SECTION  1 
INTRODUCTION 

This  final  report  has  been  prepared  by  S&D  Dynamics,  Inc.  under  contract 
DAAG29-83-C-0004  to  the  U.S.  Army  Research  Office,  Durham,  N.C.,  with  funding 
provided  by  the  U.S.  Army  Ballistic  Research  Laboratory,  Aberdeen  Proving  Ground, 
MD. 

In  accordance  with  the  contract  scope-of-work,  this  report  documents  our 
development  of  a  general,  six  degree-of-freedom  model  of  a  projectile  of  finite 
geometry  and  inertia  traveling  in  a  flexible,  rifled  or  smoothbore  gun  tube. 

The  model  accommodates  projectile  spin,  mass  eccentricity,  projectile/bore 
interfacial  friction,  elastic/plastic  deformation  of  the  projectile  rotating 
band,  and  balloting — including  bourrelet  impact  and  rebound  with  the  bore.  In 
addition,  the  model  has  been  formulated  in  terms  of  gun  tube  motion  parameters, 
thereby  rendering  the  ability  to  assess  the  mutual  effects  of  projec tile/gun- 
tube  interaction  during  in-bore  motion  when  incorporated  within,  and  solved 
simultaneously  with,  the  equations  of  a  compatible  gun  dynamics  simulation, 
such  as  DYNACODE-G  (Refs.  1,2).  Furthermore,  incorporation  of  this  model 
within  DYNACODE-G  provides  refinement  of  the  latter,  by  replacement  of  the 
traveling  point-mass  projectile  description  presently  contained  therein  with  a 
physically  realistic,  interacting  projectile  description,  as  well  as  extension 
of  the  latter  to  rifled  bore  applications.  Detailed  development  of  the  model 
is  presented  in  Section  2. 

The  basic  formulation  comprising  the  model  herein  developed  has  been 
compared  with  other  projectile  in-bore  motion  formulations  and  descriptions 
appearing  in  recent  literature.  A  summary  of  this  comparison  is  presented 
in  Section  3.  As  will  be  seen,  the  model  herein  developed  appears  to  be  the 
most  general  developed  to  date  and,  consequently,  offers  the  broadest  range 
of  applicability. 

In  view  of  the  potential  applicability  of  the  model  herein  developed,  it 
was  deemed  prudent  to  attempt  validation  of  the  model  on  its  own  merits  (prior 
to  incorporation  within  DYNACODE-G),  via  correlation  of  model  predictions  with 
experimental  in-bore  projectile  motion  data  for  firings  in  an  environment 
reasonably  free  of  gun-tube  motion.  To  achieve  this  objective,  the  U.S.  Army 
Ballistic  Research  Laboratory,  and,  in  particular,  Mr.  James  0.  Pilcner  II  and 
Dr.  James  N.  Walbert  of  the  Mechanics  and  Structures  Branch  of  the  Interior 
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Ballistics  Division,  provided  in-bore  radar  doppler  data  for  firings  of  two 
distinct  projectile  designs,  from  a  specially  designed  37mm  weapon.  Detailed 
documentation  of  the  correlation  effort  performed,  and  the  excellent  theoret¬ 
ical/experimental  agreement  demonstrated,  as  well  as  pertinent  general  observa¬ 
tions,  are  presented  in  Section  4. 

The  solution  technique  required  to  incorporate  the  model  herein  developed 
within  DYNACODE-G,  and  thereby  provide  the  desired  refinement  and  extension  of 
the  latter,  is  discussed  in  detail  in  Section  5.  In  particular,  considerations 
regarding  the  integration  scheme  and  required  step-size(s)  for  integration 
(since  the  frequency  content  in  projectile  motion  is  generally  orders  of 
magnitude  greater  than  the  frequency  content  in  basic  gun  tube  motion)  are 
given  special  emphasis. 

Finally,  conclusions  based  on  our  findings  and  recommendations  for  further 
study  are  presented  in  Section  6. 


SECTION  2 

FORMULATION  OF  PROJECTILE  EQUATIONS  OF  MOTION 

In  the  interest  of  accomodating  both  rifled  and  smoothbore  firings,  the 
formulation  herein  presented  is  developed  for  a  rifled  bore  gun  tube,  with 
relaxation  the  rifling  constraint  introduced  where  necessary  to  accomodate 
a  smt-othbore  gun  tube.  In  addition,  it  is  understood  that  the  projectile 
rotating-band  is  to  be  interpreted  as  a  bore-riding  surface  in  the  smoothbore 
application. 

The  system  of  simultaneous  differential  equations  which  describes  the 
general,  six  degree-of-freedom  motion  of  a  projectile  traveling  in  a  flexible 
gun  tube  is  formulated  first  in  terms  of  the  applied  loads  and  moments  acting 
on  the  projectile  (Section  2.1).  Based  on  practical  limitations  of  projectile 
Tti'-"'!  relative  to  the  gun  tube,  consideration  is  given  next  to  simplifying 
th/:  '<.\eral  expressions  obtained  for  the  projectile  angular  velocity  and  accel- 
e>--aion  (Section  2.2).  Finally,  the  applied  loads  and  moments  acting  on  the 
projectile  are  defined  via  modeling  of  the  projectile/gun-tube  interaction 
and  interior  ballistics  processes  (Section  2.3),  thereby  completing  the  desired 
formulation.  ■ 

It  is  noted  that  the  formulation  herein  developed  is,  within  the  framework 
of  the  assumptions  introduced,  complete  to  within  the  prescription  of  gun  tube 
motion.  It  is  further  noted  that  to  fully  account  for  the  mutual  effects  of 
projectile/gun- tube  interaction,  such  prescription  may  not  be  made  independently 
but  rather  requires  the  simultaneous  solution  of  this  formulation  with  a 
compatible  gun  dynamics  simulation,  such  as  DYNACODE-G. 

2.1  Projectile  Equations  of  Motion  in  Terms  of  Applied  Loads  and  Moments 

The  projectile  is  characterized  as  follows: 

(i)  the  projectile  consists  of  a  main  body,  rotating  band  and 
bourrelet; 

(ii)  the  projectile  main-body  behaves  as  a  rigid  body  of  finite 
geometry  and  inertia ; 

(Hi)  the  rotating  band  behaves  elastically,  with  allowance  for 
deformation  in  the  plastic  regime; 

(iv)  the  bourrelet  behaves  elastically ,  with  allowance  for  impact 
and  rebound  with  the  bore; 
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(v)  eccentricity  is  permitted  between  the  projectile  c.g.  and 
geometric  center; 

(vi)  the  projectile  is  permitted  six  degrees-of- freedom  relative 
to  the  gun  tube  (whose  motion  is  unrestricted ) ,  namely,  three 
independent  rotations  of  its  main-body  about  its  c.g.  (corre¬ 
sponding  to  projectile  pitch,  yaw  and  roll  relative  to  the  gun 
tube)  and  three  independent  translations  of  its  c.g. 

The  following  reference  frames  are  introduced  to  facilitate  tracking  projec¬ 
tile  motion: 

AAA 

(i)  S  ,  with  coordinates  (x  ,y  ,z  )  and  unit  triad  (i  ,j  ,k  ),  is 

o  o  o  o  o  o  o 

defined  as  an  intermediate  reference  frame,  fixed  neither  in 
the  projectile  nor  in  the  gun  tube,  whose  origin  translates 
with  the  projectile  along  the  gun  tube  axis  and  rotates  with 
the  projectile  about  the  instantaneous  tangent  to  the  gun 

A 

tube  axis ;  iQ  is  directed  toward  the  gun  tube  muzzle,  along 

the  instantaneous  gun  tube  axis;  at  the  initial  projectile 
/\ 

position,  j  is  directed  along  the  radius  from  the  gun  tube 

A 

centerline  to  the  projectile  c.g.;  completes  the  triad; 

(ii)  S,  with  coordinates  (x,y,z)  and  unit  triad  ( i,j,k ),  is 
defined  as  a  body-fixed  reference  frame  whose  origin  is  fixed 

A 

at  the  projectile  c.g.;  i  is  directed  toward  the  projectile 
nose,  parallel  to  its  geometric  axis;  at  the  initial  projec - 

A  * 

tile  position,  j  is  parallel  to  j  ;  k  completes  the  triad; 

^  AAA 

(Hi)  S',  with  coordinates  (x',y '  ,z' )  and  unit  triad  (i',j'.k'), 

is  defined  as  an  inertial  reference  frame  whose  origin  is 

fixed  in  space  at  the  initial  position  of  the  origin  of  Sq 

and  oriented  (for  convenience  in  subsequent  applications  of 

interest)  parallel  to  the  inertial  reference  frame  defined 

in  the  gun  dynamics  simulation  code,  DYNACODE-G. 

The  instantaneous  orientations  of  S  relative  to  S  and  S  relative  to  S'  are 

o  o 

as  depicted  in  Figures  1  and  2,  respectively.  Referring  to  these  figures  and  the 
above  definitions,  it  is  noted  that  the  kinematic  relations  between  SQ  and  S' 
incorporate  gun  tube  motions,  as  well  as  two  of  the  six  degrees-of-freedom  of  the 
projectile  relative  to  the  gun  tube,  namely,  translational  motion  along  the  gun 
tube  axis  arid  spin  (both  of  which  may  either  be  prescribed  by  interior  ballistics 
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data  or  compatibly  determined  within  the  context  of  this  formulation).  Referring 

to  Figure  1 ,  the  remaining  four  degrees-of-freedom  of  the  projectile  relative  to 

the  gun  tube,  namely,  the  two  translational  displacements  y  and  z  of  the  pro- 

A  A  «•» 

jectile  c.g.  in  the  jQ-kQ  plane  and  the  two  Euler  angles  and  $,  corresponding 
respectively  to  projectile  pitch  and  yaw,  are  formulated  in  S  relative  to  SQ. 

Referring  to  Figure  2,  the  instantaneous  position  vector  of  the  origin  of  SQ 
relative  to  S',  namely  rQ,  has  components  which  are  functions  of  time,  t,  as  well 
as  distance  along  the  gun  tube  axis,  s;  that  is 

?n(s,t)  =  x'(s,t)  i'  +  y'(s,t)  j*  +  z*(s,t)  k'  .  (1) 

Hence,  noting  that  S*  is  fixed  in  space,  and  that 


d(  )  _  d (_J_  9(  ) 

dt  ~  9t  p  9s 


(2) 


where  vp  (=  ds/dt)  denotes  the  instantaneous  velocity  of  SQ  along  the  gun  tube 
axis  (i.e.,  the  projectile  velocity  as  prescribed  by  interior  ballistics  data), 
the  instantaneous  velocity  of  Sq  relative  to  S',  namely  vq  (=  dr*Q/dt),  is  given 
by 


Vs-l)  ■  [If' 


»P  If  rfc’ ♦  »P  £’]  *  -131 


Referring  to  the  Appendix,  the  instantaneous  orientation  of  S  relative  to 


S*  is  defined  in  terms  of  the  Euler  angles  (tp0,$0,cpo)  via  the  transformation 
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where. 


SL  =  cos  ip  cos  $■ 
11  o  o 


&°  =  cos  ip  sin  $  sin  cp  -  sin  \p^  cos  <p 

12  O  O  T0  T0  O 


JL°  =  cos  ip  sin  $  cos  cp  +  sin  ip^  sin  cp„ 

1 3  T0  O  T0  O  O 


l  -  sin  \p  cos  A 

21  rO  O 


=  sin  ipQ  sin  sin  <pQ  +  cos  ipQ  cos  cpQ 


i  =  sin  \p  sin  $  cos  m  -  cos  ip  sin  cp„ 

2  3  O  O  'O  O  O 


Z  =  -  sin  A 

31  O 


*32  =  C°S  Sin  ^O 


*3  3  =  C°S  ^O  C°S  ^0 


Noting  that  i  lies  along  the  instantaneous  tangent  to  the  gun  tube  axis. 


it  follows  that 


M’  -  5°  .21'-  0°  .  i!’  =  o 

3s  '  *u  ’  3s  '  *21  ’  3s  ^31 


'«’>•>  r  \  *  -  *  -  * 
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Hence,  substituting  equations  (5)  and  (7)  into  equations  (3)  and  (4)  there 
results  the  expressions  for  the  instantaneous  translational  velocity  and  accel¬ 
eration  of  S  relative  to  S'  in  the  form 
o 


v=v  i+v  j  +  v  k 
o  xQ  o  yQ  Jo  zo  o 


where, 


v  =  v  +  £°  +  £°  3y'  +  oO  3z' 

xQ  vp  Si  at  Si  at  Si  at 


/  :£°  M’+£°  3l'  +  0° 

yQ  S2  at  *22  at  *32  at 


v  :£°  5l'  +  0°  li' 

zQ  *13  at  *23  at  *33  at 


a  =  a  i  +  a„  j  +  a  k 
°  xo  °  yo  °  zo  0 


where, 


~  -  a  A  o°  3*x*  .o  3*y’  o  3*z’ 

axQ  -  ap  +  *n  IF  *21  It^  +  *31  IF 

a  -  0°  3  x'  ^  0°  Hz'  o°  3*z'  o  r  ° 

V  '  *12  IF  +  *22  I#  +  *32  IF  +  2  Vp  Lit  cos  S  cos  \ 


3$o 

~  at2  sin  <pJ 


,2  r"o 


S  cos  9n  -  sin  <pJ 


oJ  T  'p  L3s  wo  vo  Yo  ~  Is  OAU  V 


a  -  P°  32x'  ^  0°  HZ'  2°  3*2'  o  r3^0 

aZo  -  *13  IF  +  *23  I#  +  *33  IF  -  2  Vp  [3F  C0S  S  Sln  S 

3^0  3ipo  a$ 

+  if  cos  ^3  -  vp  [35-  cos  sin  +  —  cos  <?J 


ro  T  3s 


Since  S  is  constrained  to  move  with  the  gun  tube  in  the  j  -k  plane,  it 

o  O  O  7 

has  angular  velocity,  u>  ,  relative  to  S’.  The  expression  for  a>  (written 

o  o 

relative  to  S  ),  is  given  as 


■  nrw**  vniWUl  VTfc  vm  innAn.*k  V  VlAil^Vl  Vl^APyi  v.MW  ivcv*. 


C.VUV  */>i 


o>  =  to  irt  +  w„  ,1  +  10,  k_ 
o  x  o  y  o  z  o 
o  Jo  o 


where  (referring  to  the  Appendix), 


a)  =  (p  -  ib  sin  5- 
x  Yo  ro  o 
o 


to  =  3-  cos  <p  +  4»  cos  $  sin  cp 

V  o  O  O  0  0 

Jo 


o  •-•  -  3-  sin  <p  +  Uj  cos  3  cos  <p 
z  o  Yo  o  o  To 

o 


and  ipo,  3q  and  <pQ  are  obtained  applying  equation  (2). 


Noting  that  for  a  rifled  bore  gun  tube  cpQ  is  a  function  of  s  only,  and  apply¬ 


ing  equation  (2),  there  results  the  rifling  constraint 


cp  =  t  v 
o  w  p 


where  t  (=  d<p  /ds)  denotes  the  rifling  twist;  which  is  a  known  function  pre- 
vi  o 


scribed  by  the  design  characteristics  of  the  gun  tube  rifling.  Hence,  (pQ  is  a 


known  function  of  vp  for  a  rifled  bore  gun  tube. 


Considering  a  smoothbore  gun  tube,  it  is  noted  that  investigators  generally 
assume  that  <pQ  is  identically  zero  throughout  the  in-bore  motion.  However,  as 
will  be  seen,  projectile  mass  eccentricity,  initial  projectile/bore  angular 
misalignment,  as  well  as  bourrelet/bore  impact  each  give  rise  to  non-zero  values 


of  <Pq.  Hence,  <f>o  is  herein  retained  as  an  unknown  (independent  of  vp)  for 


smoothbore  applications. 
Noting  that 


and  applying  equation  (15)  to  equation  (12),  there  results  the  angular  accelera¬ 


tion  of  S  relative  to  S’,  in  the  form 
o 


%  =  “>x  io  +  “y  jo  +  ^z  ko 
o  o  o 


where  ( from  equation  (13)), 


a)  =  9  -  sin  $  -  $  cos  $ 
x  o  o  o  o  o  o 
o 


(o  =  &  cos  cp  +  9  cos  $  sin  m  -  tL  &  sin  (L  sin  m  +  id,  m 
V  O  T0  T  o  o  ■  o  T0  o  o  To  z  To 

Jo  o 


a)  r  -  ■S-  sin  9  +  cos  &  cos  cp  -  4>  •S'  sin  $  cos  <p  -  w  cp 
z„  o  To  o  o  To  o  o  o  To  y„  o 
o  o 
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Equations  (8),  (10),  (12)  and  (16)  prescribe  the  translational  and  angular 
velocities  and  accelerations  of  Sq  relative  to  S’;  with  the  rifling  constraint, 
equation  (14),  introduced  for  a  rifled  bore  gun  tube.  It  remains  to  prescribe 
the  motion  of  the  projectile  c.g.  relative  to  S  and,  finally,  its  motion 

1 J 

relative  to  S'. 

Letting  v  and  a  respectively  denote  the  translational  velocity  and 
eg 

acceleration  of  S  relative  to  SQ,  it  follows  from  Figure  1  that 


*  /N  •  A 

.  r  1  I  m 


=  y 

j 

+  z  k 

eg 

eg 

Jo 

eg  0 

«• 

A 

•  •  A 

1 

=  y 

j 

+  z  k 

eg 

eg 

Jo 

eg  0 

(18) 


Letting  a)cg  denote  the  angular  velocity  of  S  relative  to  SQ,  and  further 
specifying  that 


(D 


eg 


=  Cl),. 


1  +  0). 


J  + 


eg 


eg 


eg 


(19) 


the  components  of  w  are  prescribed  (analogous  to  equation  (13))  as  functions 

eg 

of  the  Euler  angles  (cj),-S,cp)  of  S  relative  to  Sq  in  the  form 


•  • 

a)  =  cp  -  4»  sin  $ 
x 

eg 

♦  • 

cd  =  cos  cp  +  \J»  cos  sin  9 
ycg 


-  $  sin  9  +  41  cos  $  cos  9 


(20) 
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Noting  that  SQ  is  prescribed  to  rotate  with  S,  equation  (19)  must  be  sub¬ 


jected  to  the  constraint 


gj  „  *  i  =  0  . 
eg  o 


(21 ) 


Introducing  the  transformation  from  S  to  Sq,  namely 


AAA 

lo  =  *i.  1  *  J  +  k 

A  A  A  A 

j=J,  i  +  Jl  j  +  £  k 

J0  21  22  J  23 

•A  A  /\  /N 

kO  =  *3.  1  *  J  *  k  j 


(22) 


and  noting  that  the  direction  cosines,  are  defined  in  terms  of  the  Euler 


angles  (\p,$,cp)  analogous  to  the  definition  presented  in  equation  (6),  there 
results  from  equations  (19)  thru  (22) 


*  = 


(23) 


which  prescribes  the  dependence  between  the  Euler  angles  of  S  required  to  satisfy 
the  constraint  presented  in  equation  (21). 

Equations  (18)  and  (19)  (subject  to  equation  (23))  respectively  prescribe 


the  translational  and  angular  motions  of  the  projectile  c.g.  relative  to  SQ.  It 


remains  to  prescribe  this  motion  relative  to  S'. 

Referring  again  to  Figure  2,  the  instantaneous  position  vector  of  S  relative 
to  S’,  namely  r,  is  given  by 


r  =  r  +  r 

o  eg  * 


(24) 


Hence,  noting  equation  (15),  the  (absolute)  translational  velocity  and 
acceleration  of  the  projectile  c.g.  (relative  to  S'),  namely  v  and  a,  respec¬ 
tively,  are  given  by 


v=v  + v  +u  xr 
o  eg  o  eg 


a  =  a  +  a  +  2  ui  x  v  +  or  x  r _ +  oix  (w  x  r  ) 

o  eg  o  eg  o  eg  o  o  eg 


(25) 
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Substituting  equations  (10),  (12),  (16),  and  (18)  into  equation  (25)  and, 
in  particular,  into  the  expression  for  a  (which  is  required  for  application  to 
the  equations  of  motion  of  the  projectile) ,  there  results 


i  =  [*xo  +  2  “yo  5cg  -  2  \  ycg  *  *cg  S0  +  \  \]  -  lycg  *  el 


•  i  a  r  ••  *  •  .  • 

(a)  -o)  to  )  i  +  a  +  y  -  2  u  z  -  z _ (to 

zo  Xo  yo  J  °  L  yo  cg  xo  cg  g 


„  zcg  -  *cg  K0  -  “20  «,„> 


-  Iycg  * e)  K_  *  4.’]  3o +  [\  ♦ 2  \  hB  +  lycg  * el 


0  o 


x  (o)  +o)  a)  )  -  z  (ajz  +  a)2  )1  k 

V  f.  V  rtf  Y  “'v  J  I 


x  z  y  cg  x  y  J  o 

o  o  o  o  o 


The  (total)  angular  velocity  of  the  projectile  (relative  to  S’)  is  given  by 


0)  =  0)  +0) 

o  cg 


Applying  equation  (15)  to  each  vector  on  the  right-hand  side  of  equation 
(27)  there  results  the  (total)  angular  acceleration  of  the  projectile  (relative 
to  S') 


0)  =  0)  +  (0  +  (0  X  (0 

o  cg  o  cg 


Transforming  equation  (12)  from  SQ  to  S  via  equation  (22)  and  substituting 


the  result,  along  with  equation  (19),  into  equation  (27),  there  results 


(o=w  i+oo  j+u)  k 
x  y  z 


where , 


o)=(o  +  2,  to  +  2,  (o  +8,  (0 

x  x  n  x  21  yn  3  1  z 

cg  o  o  o 


oi=o)  +2,  0)  +2.  o)  +  2,  (o 

y  y  12  x  22  yn  32  z^ 

eg  o  o  o 


o  Jc 


uj*.  j-w-  uv-iivi-vxrv  .1.  \w\mjsajr*  ~r»  i_^.\.%r* 


Transforming  equation  (16)  from  Sq  to  S  via  equation  (22),  differentiating 
the  components  of  equation  (19)  with  respect  to  time,  and  substituting  these 
results,  along  with  equation  (19),  into  equation  (28),  chere  results 

“  =  wx  i  +  iy  j  +  k  (31 ) 


where, 


“x  =  \g  *  ln  \  ♦  h,  V  *  *„  V  -  V-  U1S  V  *  11 


eg 


w 

23  V 

O  *0 


*  *■>>  \  1  -  (».,  «,  *  H22  “y  *  I,,  W  ) 

o  eg  o  *0  o 


“y  =  “y  *  *»  \  *  l22  V  ♦  l!2  “z  *  “x  <*1,  “x  *  1 

'-B  O  O  O  Cg 


0).. 


X  '23  ~y 
O  Jo 


+  ^33  WZ  }  -  Wz  <*1,  Wx  +  *21  Wy  +  *31  CO  ) 
o  eg  o  *o  o 


“z  =  “z  *  *13  “x  -  *23  “y  ♦  *„  i,  - 


eg 


33  ~Z  “X  U12  WX  +  £22  Wy 
eg  o  yc 


*  l32  \  1  *  “y  <*„  “x  ♦  *21  1  ,  *  *  .  “x  > 
o  Jcg  o  21  o  31  zo 

and  (from  equation  (20)), 

•  ••  ••  •  • 

wx  =  cp  -  4»  sin  $  -  ij>  &  cos  & 

eg 

*  *•  ••  •  « 

0)  =  $  cos  cp  +  cos  $  sin  q>  -  $  sin  &  sin  <p  + 

Jcg 


(32) 


w2  <P 

eg 


U). 


L  =  -  $  sin  tp  +  ty  COS  $  COS  <p  -  ty  $  Sin  $  COS  m  -  (0  m 

cg  ycg 


where  (from  equation  (23)), 


(33) 


••  rtanjjj-,  V  $  r 
V  =  l  cos  $-1  $  +  ‘HST?  L~  2  ^  tan  4,  tan  $]  . 


COS  LC0S2  \|j 


(3A) 
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■-1  nr w  \  *  u  ^  w  w-i.  #.a  ^  w  u  »..*.  i  f  v  '  k "  1;  1  t.  l,  * 


Vt 


$ 


K- 


I 


1A 


r- 


R 


k 


E 

j 

Pa 


fc 

H* 


« 

K 


B-"1- 

■v 


J 


h* 
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We  are  now  in  a  position  to  formulate  the  equations  of  motion  of  the  pro¬ 
jectile  in  terms  of  the  applied  loads  and  moments. 

Letting  F  (which  remains  to  be  prescribed)  denote  the  resultant  applied 
force  acting  at  the  projectile  c.g.,  and  applying  Newton's  law  of  motion  with  a 
as  prescribed  in  equation  (26),  there  results  the  three  scalar  equations  of 


translational  motion  of  the  projectile  c.g.  (written  relative  to  SQ)  in  the  form 


F  =  m  [a  -  2  y  +  2  <d.  z  -  (y  +  e)  (u  -  u)  o>  )  ] 

v  r\  L  v  7  -  t r  /■*  rf  J  /*/*  it  v  xt  I 


p  "■  X 


z0  cg 


yQ  eg  Jcg 


x  y 
o  Jo 


V  \V] 


F  =  m  ia  +y  -2oj  z 
y0  P  y0  ycg  x0  cg 


(y„„  +  e)  (w2  +  co2  ) 
CS  x 

o  o 


-  Z  (d)  -  (0  (i)  )  1 

cg  xo  yo  zo 


F  =  m  [a  +  z  +  2  to  y„  +  (y„„  +  e)  (or,  +  w„  u)_  ) 
zo  p  2o  Cfi  xo  cg  cg  xo  yo  zo 


-  Z  ((4  +  (D2  )] 

eg  xQ  yQ 


(35) 


where  m  denotes  the  mass  of  the  projectile,  and  F  ,  F  and  F  denote  the 
p  x  y  z 

components  of  F  along  the  respective  axes  of  SQ. 


Letting  M  (which  also  remains  to  be  prescribed)  denote  the  resultant 
applied  moment  acting  about  the  projectile  c.g.,  and  applying  the  principle  of 
angular  momentum  with  u>  and  ui  prescribed  respectively  in  equations  (29)  and  (31), 
there  results  the  three  scalar  equations  of  angular  motion  of  the  projectile 
about  its  c.g.  (written  relative  to  S  in  order  to  preclude  introducing  time 
derivatives  of  the  projectile  inertia  tensor)  in  the  form 


M=I  o)+(I  -I  )  to  U)  +  I  (w  to  —  (i)  )  —  I  (u) 
x  xx  x  zz  yy  y  z  xy  x  z  y  xz  z 


+  (l)  0)  )  -  I  (d)2  -  (i)2) 

x  y  yz  y  z 
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(36) 


M=I  w  +  (I  -  X  )  w  0)  +  I  (to  to  -  to  )  -  I  (to 

y  yy  y  xx  zz  x  z  yz  x  y  z  yx  x 

+  to  to  )  +  I  (to2  -  w2) 
y  z  zx  x  z' 

•  •  • 

M  =  I  to  -  (I  -  I  )  to  to  +  I  (w  to  -  to  )  -  I  (to 

z  zz  z  xx  yy  x  y  zx  y  z  x  zy  y 

+  w  w  )  -  I  (to2  -  to2) 
x  z  xy  x  y 

where  M  ,  M  and  M  denote  the  components  of  M  about  the  respective  axes  of  S, 
x  y  z 

and  I  ,  I  I  .  denote  the  elements  of  the  projectile  inertia  tensor 

xx’  yy’  zy 

written  relative  to  S. 


2.2  Simplification  of  Angular  Velocity  and  Acceleration  Expressions 

The  expressions  for  the  terms  representing  the  angular  velocity  and  accel¬ 
eration  components  entering  equation  (36)  may  be  greatly  simplified  by  noting 
that  fcr  most  practical  applications  \|>  and  0  are  sufficiently  small  such  that 


sin  9  a  9  ;  cos  ~  1 


sin  $  ss  $  ;  cos  $  «  1 


Under  this  condition  it  follows  from  equation  (23)  that 


9  *  ^  &  . 


1 


Applying  the  mean  value  theorem,  it  follows  that  the  integral  of  <p  will  at 
most  be  of  order  9  -0 .  Hence,  to  first  order 


cp  ~  0 


Imposing  the  above  conditions  and  retaining  only  linear  (first  order) 
terms  in  \}>  and  3,  the  direction  cosines  for  the  transformation  from  S  to  SQ, 
as  given  in  equation  (22),  simplify  to 


£11  =  £22  =  £33  =  1 


£  =  -  £  =  -  U> 

12  21 


£13  =  -*31  »* 


£23  =  £32  =  0 


Noting  the  above,  equations  (20)  and  (33)  simplify  respectively  to 


(o  =  i|j  -9-  —  -9-  vp 
Xcg 


U,  *9 

eg 


WZ  *  * 

eg 


a)  ~  tjj  9  -  9  i}; 

Xcg 


a)  »  9  +  i)>  ip  9 
Ycg 


to  »ip-994>-4)92 
z 

eg 


Substituting  equations  (40)  thru  (42)  into  equations  (30)  and  (32),  there 
results  the  greatly  simplified  expressions  for  the  angular  velocity  and  accel¬ 
eration  components  entering  equation  (36)  in  the  fora 


wx  =  wx  +  ^  ^wv  +  ^  ^  +  ^ 

xo  yo  zo 


wv  =  w  +  9  -  4>w 
y  yo  o 


wz  =  “z  +  ♦  +  *  Wx 
o  o 


and 


a)  =  a)  +  (^  +  w  )  -  $  ("i|>  +  o>  )  +  $  ($  ul,  +  w _  ) 

xx  y  z  x  z 

O  J0  0  o  o 

+  4>  (\j)  0)  -  Cd  ) 

o  yo 

to  =  to  +  $  +  i|3  (•&■’!’-<*>  +  t|3  w  ) 

y  yo  xo  xo  yo 

+  il>  &  {$  w  +  w  ) 

X  z 
o  o 

•  •••  •  •••  • 

03  =  0)  +  l}/-\J/'9-2  -  •&{$•  l|)-U3  )+  &  (0)  -  &  0)  ) 

2  X  X  Z 

o  o  o  o 

-  ■&  vp  (ib  o)  -  0)  ) 

x  y 

A  •'A 


(44) 


2.3  Applied  Loads  and  Moments 

The  equations  of  projectile  motion  developed  in  the  preceding  section 
require  specification  of  the  applied  loads,  F  (written  relative  to  Sq)  ,  and 
the  applied  moments,  M  (written  relative  to  S),  for  completion  of  the  formula¬ 
tion.  These  loads  and  moments  arise  as  a  consequence  of  the  projectile  weight, 
interfacial  contact  of  the  rotating  band  and  bourrelet  with  the  bore,  propel¬ 
lant  gas  pressure  acting  at  the  base  of  the  projectile,  and  compressed  air 
ahead  of  the  projectile. 


2.3.1  Projectile  Weight  Loading 

The  load  applied  to  the  projectile  c.g.  due  to  its  own  weight  is 
given  by 

F  =  -  mn  g  3f  (*5) 

W  p 

where  m  and  g  denote  respectively  the  projectile  mass  and  gravitational 

p 

acceleration. 

Applying  the  transformation  given  in  equation  (5),  equation  (45)  is 
written  relative  to  SQ  for  application  to  equation  (35)  in  the  form 
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F  _  F  i  +  F  1  +  F  k 

w  "  xn  o  v  Jo  o 

O  O  0 


tn  g  sin  cos  $ 
p  o  o 


m  g  (sin  ^  sin  sin  cp  +  cos  cos  cp  ) 
p  o  o  o  o  o 


m  g  (sin  i>  sin  $  cos  cp  -  cos  4<  sin  cp  ) 
p  O  O  0  o  o 


2.3.2  Rotating-Band/Bore  Interfacial  Contact  Loading 


The  loads  and  moments  applied  at  the  projectile  c.g.  due  to  interfacial 
contact  of  the  rotating  band  with  the  bore  are  derived  subject  to  the  following 
assumptions: 

(i)  since  modeling  the  engraving  process  in  a  rifled  bore  gun  tube 
is  beyond  the  scope  of  this  effort,  the  initial  state  of  the 
projectile  and,  in  particular,  the  rotating  band  is  defined 
here  as  its  fully  engraved  state;  hence,  for  this  purpose, 
the  initial  state  of  the  rotating  band  is  characterized  as 
a  radially  compressed  elastic  band  with  elastic/plastic 


boundary  at  the  bore  surface ; 


(ii)  projectile  motion  subsequent  to  engraving  induces  additional 
elastic,  as  well  as  plastic  radial  deformations  of  the 
rotating  band; 

(Hi)  elastic  radial  deformation  of  the  rotating  band  is 

characterized  by  a  Winkler  foundation  model;  rendering 
a  radially  directed  rotating-band/bore  interfacial  load 
distribution  (with  local  magnitude  determined  by  resultant 
local  radial  displacement)  which  varies  circumferentially 
around  the  bore,  as  well  as  longitudinally  along  the 
length  of  the  interface;  the  latter  variation  gives  rise 


to  the  so  called  " foundation "  moment; 


( iv)  each  cross-section  of  the  rotating  band  maintains  full 


contact  with  the  bore; 


(v)  Coulomb  friction  acts  at  the  rotating-band/bore 
interface;  rendering  a  load  distribution  with  magnitude 
proportional  to  the  radially  directed  distribution 

and  directed  opposite  to  the  resultant  motion  at 
the  interface; 

(vi)  rifling  torque  is  transmitted  without  slippage  via  a 
uniformly  distributed  circumferential  load  acting  along 
the  length  of  the  rotating  band  at  the  rotating-band/ 
bore  interface; 

(vii)  the  projectile  main-body  is  rigid  compared  to  the 
rotating  band; 

(viii)  the  Euler  angles  ^  and  -8-  satisfy  equation  (37). 

In  view  of  Assumptions  (vii)  and  (viii),  the  displacement  relative  to 
SQ,  namely  6,  of  any  point  of  the  projectile  in  a  plane  perpendicular  to  its 
geometric  axis  is,  at  any  instant  during  the  motion,  given  by 


« *  <ycg  -  5  *>  J0 *  <zcg <«> 

where  £  denotes  the  perpendicular  distance  from  the  projectile  c.g.  to  the 
plane  of  interest.  It  is  noted  that  £  >  0  implies  that  the  plane  of  interest 
is  aft  of  the  projectile  c.g.;  £  =  0  implies  that  the  plane  of  interest  con¬ 
tains  the  projectile  c.g.;  while  £  <  0  implies  that  the  plane  of  interest 
is  forward  of  the  projectile  c.g. 

The  radial  component  of  this  displacement,  namely  6r,  is  given  as 

—  A 

6  =  6  •  i  =  (y  -  £  \Jj)  cos  <p  +  (z  +  £  $)  sin  <p  (49) 

I  i  eg  eg 

where  cp  denotes  the  angle  between  the  projection  of  the  y  -axis  onto  the  plane 
of  interest  and  the  line  from  the  gun  tube  centerline  to  the  point  under 
consideration. 

The  maximum  radial  displacement,  6  ,  in  the  plane  of  interest  is 

max 

determined  by  setting 
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V- v> 

iV-A 


.  -.w 


6=6  cos  (cp  -  cp  )  (53) 

which  prescribes  the  circumferential  distribution  of  radial  displacement  in 

the  plane  of  interest.  It  is  noted  that  this  distribution  is  symmetric  with 

respect  to  the  angle  cp  .  It  is  further  noted  that  both  6  and  cp  vary 

iDuX  max  max 

from  one  plane  to  the  next  in  accordance  with  the  variation  of  £.  Hence, 
applying  equations  (51)  thru  (53)  to  planes  containing  the  rotating  band,  it 
is  seen  that  the  radial  deformation  within  a  band  cross-section  is  non-uni- 
formly  distributed  around  its  circumference,  and  that  both  the  magnitude  and 
orientation  of  this  distribution  vary  from  one  cross-section  to  the  next  along 
the  length  of  the  band. 

Superposing  the  initial  radial  compression,  6Q,  in  accordance  with 
Assumption  (i),  and  multiplying  by  the  elastic  ’’spring"  stiffness  per  unit 
surface  area  of  the  rotating  band,  k,  there  results  in  accordance  with 
Assumption  (iii),  the  radial  load  per  unit  surface  area,  R(cp,£),  acting  on 
the  rotating  band  at  the  rotating-band/bore  interface  at  each  instant  during 
the  motion,  in  the  form 

R(«P,5)  =  -  k  (6p  +  6q)  ip  •  (54) 

In  view  of  Assumption  (iv),  the  radial  load  distribution  at  a  typical 
cross-section  of  the  rotating  band,  as  specified  by  equation  (54)  for  a 
constant  value  of  £,  is  as  depicted  in  Figure  3. 


In  view  of  Assumptions  (iv)  thru  (vi),  the  instantaneous  interfacial 
contact  loads  acting  on  a  differential  element  of  the  rotating  band  are  as 
depicted  in  Figure  4. 

Under  Assumption  (v),  the  friction  force,  pi?,  acts  tangential  to  the 
bore  surface  in  the  direction  opposing  projectile  motion  at  the  surface.  The 
angle  a  depicted  in  Figure  4  is  determined  from  the  expression 


tan  a  =  < 


r.  t  ;  rifled  bore 
b  w  ’ 


rb  *0 
v  VP 


;  smoothbore 


(55) 


In  addition,  in  accordance  with  Assumption  (vi),  the  rifling- torque  load 
per  unit  surface  area  of  the  rotating  band,  T,  is  constant  over  the  entire 
rotating-band/bore  interface  in  a  rifled  bore  gun  tube;  whereas,  T  =  0  for  a 
smoothbore  gun  tube. 

Letting  dF  denote  the  resultant  incremental  contact  load  acting  on  a 

differential  element  of  the  rotating  band,  it  may  be  written  in  component  form 

(relative  to  S  )  as 
o 


dF. 


yQ  J°  * 


dF 


(56) 


where  referring  to  Figure  4 

dF^  =  -  rfa  p  R(<p,5)  cos  a  d£  dcp 
o  u 

<p]  -  T  sin  (p}  d£  dtp 

<p]  +  T  cos  cp}  d?  dtp 
o 


dFc  =  rfa  {R(<p,  S)  [-  cos  <p  +  p  sin  a  sin 
^o 

dFc  =  r.  {r  { tp,  £)  [-  sin  tp  -  p  sin  a  cos 

Z  D 


)  (57) 


Integrating  over  the  rotating-band/bore  interfacial  surface  area,  there 
results,  for  application  to  equation  (35),  the  contribution  to  the  force 
components  acting  at  the  projectile  c.g.  (relative  to  Sq)  due  to  rotating 
band  contact  with  the  bore,  as  follows 
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(58) 


F  :  *  2  n  p.  p  k*  i  cos  a 
x  bo 

o 

Fy  =  "  rb  k'  l(lb  '  I1  *  "  ycg  *  w  si"  “  tzcg  *  (ib  ■  I’ 

Fz  =  -  "  rb  k'  lzcg  ♦  ab  -  t>  4  ♦  >*  sin  a  [Fcg  -  (1b  -  f>  «) 


where  denotes  the  distance  from  the  projectile  c.g.  to  the  rear  face  of  the 
rotating  band  (which  is  assumed  to  be  located  aft  of  the  c.g.),  h  denotes  the 
width  of  the  rotating  band,  and  k'  (=  k  h)  denotes  the  "spring"  stiffness  of 
the  rotating  band  per  unit  circumferential  length. 

Referring  again  to  Figure  4,  the  moment  arm,  from  the  projectile  c.g. 
to  the  load  acting  on  the  differential  element  is  given  by 

V  -  U0  +  [rb  cos  <p  -  (ycg  +  e>]  J0  +  [rb  sin  <p  -  zQg]  kQ  .  (59) 

Hence,  the  resultant  moment  applied  to  the  projectile  c.g.  is  given  by 

M  =  /;  L  x  d?  (60) 

where  denotes  the  instantaneous  rotating-band/bore  interfacial  surface  area. 

Substituting  equations  (56),  (57)  and  (59/  into  equation  (60),  and  perform¬ 
ing  the  indicated  integration,  there  results  the  contribution  to  the  moment 
components  acting  at  the  projectile  c.g.  (relative  to  SQ)  due  to  rotating  band 
contact  with  the  bore,  as  follows 

Mx  *  rb  "  rb  »  T  -  tan  “  F^  5  *  zcg  FJ  -  (ycg  *  E>  Fz 

o  o  o  o 


M  =  -  z 


Ff  +  (A.  -  §)  F®  -  ir  r2  k»  V  cos  a  (z 


eg  x  b  2'  z 

o  o 


7T  v  h2  k  * 

+  (£b  -  |)  + - % - [  (y  sin  a)  ij)  -  $] 


MC  =  (y  +  c)  FC  -  (£.  -  FC  +  tt  r2  k*  M  cos  a  [y 
zq  Jcg  xq  b  2  yQ  b  LJcg 

tt  r  h2  k* 

-  <*b  -  |)  *] - J2 -  £*  +  (p  sin  a)  ^ 


i  .2  - 


where  Fc  ,  Fc  and  FC  are  as  defined  in  equation  ( 58 ) ,  and  T  =  0  for  a 
xo  yo  zo 
smoothbore  application. 

Referring  to  equation  (61),  it  is  of  interest  to  note  that  characterization 

of  the  elastic  response  of  the  rotating  band  by  a  Winkler  foundation  model  gives 

rise  to  the  so  called  "foundation"  moment,  which  manifests  itself  via  the  last 

c  c 

term  on  the  right-hand  sides  of  the  expressions  for  M  and  M  . 

yo  zo 

The  desired  contribution  to  the  moment  components  M  ,  M  and  M  written 

x  y  z 

relative  to  S  for  application  to  equation  (36),  is  obtained  applying  the  inverse 
of  the  transformation  given  in  equation  (22),  from  which  there  results  (noting 
Assumption  (viii)) 


MC  =  MC  +  MC  -  $  MC 
x  xo  yo  Z, 

MC  =  MC  -  t{)  MC 
y  yo  Xo 

MC  =  MC  +  $  MC 

Z  Z  X 


(62) 


2.3.3  Bourrelet/Bore  Interfacial  Contact  Loading 

To  accommodate  a  variety  of  bourrelet  designs,  we  distinguish  between 
bourrelets  which  are  initially  either  bore  diameter  or  greater,  or  sub-caliber. 
In  either  case,  to  facilitate  computations  allowing  for  bourrelet/bore  contact, 
the  instantaneous  curvature  of  the  gun  tube  axis  between  the  planes  containing 
the  rear  face  of  the  rotating  band  and  the  forward  face  of  the  bourrelet  is 
neglected. 

Considering  the  case  wherein  the  bourrelet  is  either  bore  diameter  or 
greater,  it  is  treated  as  was  the  rotating  band  in  the  preceding  section. 

Hence,  the  equations  developed  in  the  preceding  section  and,  in  particular, 
equations  (58),  (61)  and  (62),  with  replaced  by  the  negative  of  the  distance 
from  the  projectile  c.g.  to  the  forward  face  of  the  bourrelet  (assuming  the 
bourrelet  is  forward  of  the  projectile  c.g.),  h  replaced  by  the  width  of  the 
bourrelet,  and  y  and  k*  replaced  by  appropriate  values  for  the  bourrelet  remain 
valid  for  this  case. 

Considering  the  case  wherein  the  bourrelet  is  sub-caliber,  we  accommodate 
bourrelet  impact  with  the  bore  by  applying  the  principles  of  linear  and  angular 
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momentum  (in  integrated  form)  to  determine  the  impulsive  motions  imparted  to 
the  projectile  during  impact.  It  is  noted  that  at  the  instant  of  impact  an 

a 

impulsive  load,  F,  is  generated  at  the  bourrelet/bore  interface  and,  in  addi¬ 
tion,  since  projectile  spin  is  constrained  by  the  rifling  twist  in  the  rifled 

A 

bore  application,  an  impulsive  torque,  T,  is  simultaneously  generated  at  the 
rotating-band/bore  interface. 

Letting  A(  )  denote  the  "jump”  (sudden  increment)  in  the  parameter  (  ) 
due  to  impact,  the  principle  of  linear  momentum  (integrated  with  respect  to 
time)  may  be  written  in  the  form 

m  Av'  =  f  (63) 

P  eg 


where  v'  denotes  the  velocity  of  the  projectile  c.g.  relative  to  the  gun  tube, 
eg 

Similarly,  the  principle  of  angular  momentum  (integrated  with  respect  to 
time)  may  be  written  in  the  form 

AH  =  rtx  F  +  f  (64) 


where  H  denotes  the  angular  momentum  of  the  projectile  about  its  c.g.  and  r 
denotes  the  position  vector  from  the  projectile  c.g.  to  the  point  of  applica¬ 
tion  of  the  impulsive  load. 

Referring  to  Figure  5  and  the  notation  introduced  in  Section  2.1 


F 


=  -  [pj  cos  a]  R  iQ  -  [cos  q>j  -  pt  sin  a  sin  cpj]  R  jQ 

A  /* 

-  [sin  cpj  +  px  sin  a  cos  R  kQ 
f  -  T  iQ  ;  rifled  bore 


>  (65) 


T  = 


smoothbore 


f>  *  5  ^  +  [rb  cos  *i  -  ,ycg  +  C)1  *  [rb  sln  *1  -  zcg]  “o 


where  R  denotes  the  radially  directed  impulsive  load  acting  at  the  bourrelet/ 

A 

bore  interface,  T  denotes  the  circumferentially  directed  impulsive  torque  acting 


v- 

iS 
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l\-\4  > 


U-i. 


Figure  5  -  Bourrelet/Bore  Impact  Loadin, 


at  the  rotating-band/bore  interface,  y  denotes  the  friction  coefficient  at  the 
bourrelet/bore  interface,  a  is  as  defined  in  equation  (55),  and  and  <pt  are 
as  depicted. 

Referring  again  to  the  notation  of  Section  2.1  and,  in  addition,  assuming 
the  validity  of  equation  (37),  it  may  be  shown  that 


v1  =  v  i  +  [y  -  z  <p  ]  j  +  [z  +(y  +  e)  m  ]  k 

eg  p  o  Jcg  eg  Yo  Jo  eg  Jcg  ro  o 

H  =  I  (cp  +  \p$--&ip)i  +  I  (■&  -  ip  cp  )  j  +  I  (\p  +  •&  <p  )  k 
xx  To  yy  Y  ro  J  zz  To 


}  (66) 


Substituting  equations  (65)  and  (66)  into  equations  (63)  and  (64),  per¬ 
forming  the  indicated  vector  operation  on  the  right-hand  side  of  equation  (64) 
and,  furthermore,  transforming  the  resulting  expressions  from  SQ  to  S  (noting 
equation  (37)),  there  results  the  six  scalar  equations 


mp  Avp  =  -  (y1  cos  a)  R 

m  (Ay _  -  z _  Am)  =  (-  cos  m  +  y,  sin  a  sin  o  )  R 

p  eg  eg  To  TI  1  Ti 

m  [Az  +  (y  +  e)  Aq>  ]  =  -  (sin  <j>  +  y  sin  a  cos  <p  )  R 
p  eg  eg  O  11  1 

Ixx  +  ip  AS  -SAip)  =  R  {e  sin  tpi  +  yt  sin  a  [(ycg  +  e 

-  ip)  cos  <pt  +  (zcg  -  &i  S)  sin  <pi  -  rfa]  -  yt  cos  a 

x  [rb  (ip  sin  cpi  +  $  cos  <pt )  -  e  S]}  -  T 

I  (AS  -  ip  A(pQ)  =  R  {(g>i  -  e  vp)  sin  cpt  +  yt  sin  a  [(2^  -  e  ip) 

1  A 

x  cos  <p.  +  r.  ip]  -  y  cos  ot  (r.  sin  cp  -  z  )j  +  \p  T 
Ti  d  i  b  Ti  eg 

Izz  ^  +  ^  A<pQ)  =  R  1“  cos  <Pi  +  e  9  sin  +  Vij  sin  a 


>  (67) 


x  (il  sin  <pj  +  e  S  cos  cpj  -  rfa  S)  +  yj  cos  a  [rfa  cos  cpj 
-  (ycg  +  e)3}  -  3  T 
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i 

A 

with  T  =  0  in  the  smoothbore  application. 

In  addition,  the  constraint  imposed  by  the  rifling,  namely  equation  (14), 
requires  that 

to  :  t  4v  (68) 

To  w  p  • 

Hence,  considering  the  rifled  bore  application,  equations  (67)  and  (68) 
represent  seven  algebraic  equations  in  the  eight  unknowns,  Av  ,  Ay  ,  Az  , 

•  •  •  A  A  POO 

AcpQ,  A$,  A4>,  R,  and  T;  whereas,  considering  the  smoothbore  application,  equa¬ 
tion  (67)  contains  seven  unknowns.  The  remaining  equation  required  to  complete 
the  set  for  either  application  is  obtained  by  characterizing  the  impact  process 
at  the  bourrelet/bore  interface  via  introduction  of  the  coefficient  of  restitu¬ 
tion,  e. 

In  accordance  with  the  above  notation  and  assumptions,  the  velocity  rela¬ 
tive  to  the  bore  of  the  point  on  the  bourrelet  which  impacts  the  bore  is  given 
by 

_  .  ..  A 

v’  =  {v  +  (r.  sin  <p  -  z  )  &  -  [r.  cos  cp  -  (y _ +  e)3  ip}  i 

c  1  p  b  ri  eg  b  Yj  eg  TJ  o 

+  +  a  -  rK  <p~  sin  <p  1  +  iz^„  -  i  & 

eg  i  Y  b  Yo  YiJ  o  1  eg  i 

•  ^ 

+  r.  cn  cos  cp  }  k  (69) 

b  Yo  Yi 1  o  • 

Since  compression  and  restitution  during  impact  occur  in  a  direction  normal 

A  _ 

to  the  bore  surface  (along  the  line  of  action  of  R),  the  component  of  v'  normal 
to  the  bore  surface  immediately  following  restitution,  namely  (v^»  i )  ,  is 
related  to  its  value  just  prior  to  impact,  namely  (v£  •  ip)~,  via  the  coeffi¬ 
cient  of  restitution,  that  is 


( v'  •  i)+  =  -  e  (v'  •  i  )"  .  (70) 

c  r  c  r 

Hence,  there  results  from  equations  (69)  and  (70) 

(Aycg  +  Acj))  cos  <pi  +  (Azcg  -  di  AS)  sin  cpt  =  (1  +  e) 

*  _  •  —  _  •  _ 

x  [(y  +  l  )  cos  cp  +  (z  -  fi.  $  )  sin  cp  ]  .  (71) 

eg  i  Y  Yi  eg  1  Yi 
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Equations  (67),  (68)  and  (71)  in  the  rifled  bore  application,  and  equations 
(67)  and  (71)  in  the  smoothbore  application,  completely  prescribe  the  "jump" 
conditions  at  impact  in  terms  of  projectile  orientation  and  motion  just  prior 
to  impact  and  the  relative  material  properties  of  the  bourrelet  (via  the 
coefficient  of  restitution).  It  is  noted  once  again  that,  in  general,  bourrelet 
impact  with  the  bore  gives  rise  to  non-zero  increments  in  the  parameter  for 
both  rifled  and  smoothbore  applications. 

2. 3. A  Propellant  Gas  Pressure  and  Compressed  Air  Loadings 
The  loads  and  moments  applied  to  the  projectile  c.g.  due  to  propellant  gas 
pressure  acting  at  the  base  of  the  projectile  and  "ram"  (compressed)  air  ahead 
of  the  projectile  are  derived  subject  to  the  following  assumptions: 

(i)  both  the  projectile  base  and  "ram"  air  pressures  are  assumed  to 
be  known  functions  of  time  only,  and  at  any  instant  are  uniformly 
distributed  over  the  respective  projectile  surfaces  over  which 
they  act; 

(ii)  for  cases  wherein  the  bourrelet  is  either  bore  diameter  or 
greater,  both  the  rotating  band  and  bourrelet  are  assumed  to 
act  as  ideal  gas  seals ;  for  cases  wherein  the  bourrelet  is 
sub-caliber,  the  rotating  band  alone  is  assumed  to  act  as  an 
ideal  gas  seal; 

(Hi)  the  Euler  angles  i|>  and  $  satisfy  equation  (37). 

It  follows  from  Assumption  (i)  that  the  base  and  "ram"  air  pressures  act 

/\ 

as  effective  hydrostatic  pressures;  generating  a  force  equal  to  -  (p  dA)  i 
acting  on  a  differential  surface  area  element,  dA,  with  unit  outward  normal 

A 

in«  Further,  from  general  principles  of  hydrostatics  and  Assumption  (ii), 

these  pressure  loadings  may  be  replaced  by  resultant  forces  acting  at  the 

geometric  center  of,  and  directed  perpendicular  to,  the  sealing  planes. 

Considering  first  cases  wherein  the  bourrelet  is  either  bore  diameter  or 

greater,  the  sealing  planes  contain  the  forward  face  of  the  bourrelet  and  the 

rear  face  of  the  rotating  band.  Hence,  the  points  of  application  of  the 

resultant  pressure  loads  are  defined  by  the  respective  intersection  of  each 

of  these  planes  with  the  gun  tube  centerline,  as  depicted  in  Figure  6. 

/\ 

Referring  to  Figure  6,  and  noting  that  since  i  is  parallel  to  the  geometric 
axis  of  the  projectile,  it  follows  that 
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Figure  6  -  Pro.jectile  Base  and  "Ram"  Air  Pressure  Loads 


where  p  and  p  denote  respectively  the  instantaneous  "ram"  air  and  projectile 

a  D 

base  pressures,  and  A  denotes  the  projected  area  of  the  bourrelet  and  rotating 

A 

band  planes  perpendicular  to  i. 

Under  Assumption  (iii),  it  follows  that 


A  *  ir  r£ 


(73) 


Hence,  the  resultant  pressure  load,  Fp,  is  given  by 


Fp  *  *  rb  lpb  -  pa’  1 


(74) 


which  is  written  relative  to  SQ  for  application  to  equation  (35)  in  the  form 


F  =  fP  i  +  FP  i  +  FP  k 
p  xo  °  \  °  \  ° 


(75) 


where  (noting  Assumption  (iii)), 


Fx  =  ’  rb  ,pb  -  pa! 
o 


Fyo  =  "  rb  (pb  *  pa'  * 


Fz  =  ”rb  (pb  -  pa’  9 
o 


(76) 


Once  again  in  view  of  Assumption  (iii),  the  moment  arm,  SL  ,  from  the 

Si 

projectile  c.g.  to  the  "ram"  air  load,  F  ,  is  given  by 

Si 


\  =  ft,  i  -  (ycg  +  e  +  1|| )  j  -  (zcK  -  ai  $)  k 


eg 


(77) 


where  l1  denotes  the  distance  from  the  projectile  c.g.  to  the  forward  face  of 
the  bourrelet. 
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The  corresponding  moment  arm,  £b,  to  the  load  Fb  is  given  by 


*b  =  •  *b  1  -  lycg  *  e  -  *b  ♦'  J  -  l2cg  *  *b  41  k 


where  &b>  as  previously  defined,  denotes  the  distance  from  the  projectile  c.g. 


to  the  rear  face  of  the  rotating  band. 

Hence,  the  resultant  moment  applied  to  the  projectile  c.g.  is  given  by 


M  =  S,  x  F  +  2,.x  F.  . 
p  a  a  b  b 


Substituting  equations  (72),  (77)  and  (78)  into  equation  (79),  and  per¬ 
forming  the  indicated  vector  operations,  there  results  the  desired  contribu¬ 
tion  to  the  moment  components  acting  at  the  projectile  c.g.  (relative  to  S) 
due  to  propellant  gas  and  "ram”  air  pressures,  in  the  form 


MP  =  0 
x 


MP  =  -  n  r£  [(p  -  p  )  z  +  &  (p.  i .  +  p  l  )] 

y  b  b  a  eg  b  b  a  i 


Mz  2  *  rb  t(pb  -  pa’  (ycg  *  El  -  *  lpb  lb  *  pa  V] 


Considering  the  case  wherein  the  bourrelet  is  sub-caliber,  equation  (76) 
remains  unchanged,  while  equation  (80)  is  modified  by  replacing  by 


(Jlb  -  h) ;  where  h,  as  previously  defined,  denotes  the  width  of  the 


rotating  band. 


2. A  Solution  Technique 

Summing  corresponding  force  components  from  equations  (47),  (58)  and  (76), 
there  results  for  application  to  equation  (35) 


F  -  fW  +  FC  +  FP 

X  X  X  X 

O  O  O  0 


F  =  FW  +  FC  +  FP 

y  y  y  y 

Jo  Jo  Jo  Jo 


F  -  FW  +  FC  +  FP 
z  z  z  z 
O  0  o  o 


Similarly,  summing  corresponding  moment  components  from  equations  (62)  and 
(80),  there  results  for  application  to  equation  (36) 


M  =  MC  +  MP 


M  =  MC  +  MP 
z  z  z 


(82) 


Within  the  framework  of  the  assumptions  introduced,  equations  (14),  (35) 
and  (36),  with  the  applied  loads  and  moments  as  defined  in  equations  (81)  and 
(82),  and  the  projectile  "jump"  conditions  due  to  bourrelet/bore  impact  as 
defined  in  equations  (67),  (68)  and  (71),  constitute  the  desired  formulation 
for  a  rifled  bore  gun  tube;  whereas,  deleting  equations  (14)  and  (68),  there 
results  the  desired  formulation  for  a  smoothbore  gun  tube. 

Considering  either  formulation,  it  is  seen  that  projectile  in-bore  motion 
is  prescribed  in  terms  of  projectile  design  data,  interior  ballistics  data,  and 
gun  tube  design  and  motion  data.  Of  these  data,  gun  tube  motion  is  generally 
not  known  a  priori.  In  fact,  an  objective  of  this  formulation  is  to  define  the 
mutual  effects  of  projectile/gun-tube  interaction.  Hence,  to  satisfy  this 
objective,  the  formulation  herein  developed  (for  both  rifled  and  smoothbore 
firings)  is  to  be  solved  simultaneously  with  a  compatible  gun  dynamics  simula¬ 
tion,  such  as  DYNACODE-G.  In  addition,  since  the  onset  of  bourrelet/bore  impact 
is  also  generally  not  known  a  priori,  the  solution  technique  requires  monitoring 
projt  le  in-bore  motion  at  each  instant  to  determine  such  onset  and,  accord¬ 
ingly  ntroduction  of  the  projectile  "jump"  conditions. 


SECTION  3 

COMPARISON  WITH  OTHER  PROJECTILE  MOTION  FORMULATIONS 

The  formulation  developed  in  Section  2  is  general  in  that  it  prescribes 
the  full  six  degree-of- freedom  motion  of  a  projectile  of  finite  geometry  and 
inertia  traveling  in  a  flexible  (rifled  or  smoothbore)  gun  tube.  The  degrees- 
of-freedom  selected  correspond  to  three  independent  translational  motions  of 
the  projectile  c.g.  relative  to  the  gun  tube  axis  and  three  independent  rota¬ 
tions  of  the  projectile  about  its  c.g.  (corresponding  to  pitch,  yaw  and  roll 
motions).  In  addition,  the  formulation  accommodates  bourrelet  impact  and 
rebound  with  the  bore  (via  impulsive  loadings  which  give  rise  to  "jump" 
conditions  in  projectile  motion  parameters  during  impact).  Finally,  since 
gun  tube  motion  (which  is  unrestricted  in  this  formulation)  is  generally  not 
known  a  priori,  the  projectile  equations  of  motion  have  been  formulated  with 
the  intent  of  being  solved  simultaneously  with  the  equations  of  a  compatible 
gun  dynamics  simulation,  such  as  DYNACODE-G. 

A  similar,  but  far  more  restrictive  formulation  has  recently  been  pre¬ 
sented  by  S.H.  Chu  (Ref.  3).  Chu  permits  the  projectile  three  degrees-of- 
freedom;  two  orthogonal  translational  motions  of  the  projectile  c.g.  relative 
to  the  gun  tube  and  one  rotational  (pitching)  motion.  As  a  consequence  of 
neglecting  the  remaining  degrees-of-freedom,  the  gun  tube  centerline,  the 
projectile  c.g.,  and  the  resultants  of  the  rotating  band  and  bourrelet  contact 
load  distributions  with  the  bore  all  lie  in  the  same  plane.  Hence,  Chu's 
formulation  is  essentially  planar;  whereas,  the  formulation  herein  presented 
is  of  general  three-dimensional  character. 

Another  recent,  but  also  restrictive  projectile  motion  formulation  has 
been  presented  by  H.L.  Langhaar  and  A.P.  Boresi  (Ref.  A).  Langhaar  and  Boresi 
present  a  rigorous  kinematical  description  of  a  point  moving  along  a  time-depend¬ 
ent  space  curve.  The  point  is  identified  with  the  geometric  center  of  a  rigid 
projectile;  the  time-dependent  space  curve  is  identified  with  the  centerline  of 
a  flexible  gun  tube.  The  rigid  projectile  is  further  characterized  such  that 
its  geometric  center  and  c.g.  coincide  (which  precludes  the  ability  to  investi¬ 
gate  the  effects  of  mass  eccentricity),  and  such  that  its  geometric  axis  is 
directed  along  the  instantaneous  tangent  to  the  gun  tube  centerline  (which 
precludes  the  ability  to  investigate  the  effects  of  projectile  pitch  and  yaw 
motions  relative  to  the  gun  tube).  The  projectile  is  permitted  two  degrees- 


42 


of-freedom  relative  to  the  gun  tube;  translational  motion  of  its  c.g.  along  the 
gun  tube  centerline  and  rotational  motion  about  the  centerline  ( correspondi^ 
to  projectile  spin).  Accounting  for  rotary  inertia  of  the  projectile  about  its 
spin  axis,  there  results  a  traveling  point-mass  projectile  load  with  superposed 
gyroscopic  couple. 

J.J.  Wu  (Ref.  5)  also  adopts  a  traveling  point-mass  projectile  description, 
but  with  superposed  traveling  pitching  moment  due  to  mass  eccentricity  (while 
neglecting  rotary  inertia  about  the  pitch  axis).  Several  other  investigators 
(Refs.  6  thru  8)  have  adopted  the  simpler  traveling  point-mass  description, 
with  and  without  mass  eccentricity  and  projectile  spin  (while  neglecting 
rotary  inertia  about  the  spin  axis). 


SECTION  4 

CORRELATION  WITH  EXPERIMENTAL  DATA 

Before  incorporating  the  formulation  herein  developed  within  DYNACODE-G, 
it  was  deemed  prudent  to  verify  its  validity  on  its  own  merits.  However,  to 
achieve  such  verification  required  establishing  an  experimental  data  base,  as 
well  as  programming  the  model  equations. 

To  establish  the  experimental  data  base,  Mr.  James  0.  Pilcher  II  and 
Dr.  James  N.  Walbert  of  the  Mechanics  and  Structures  Branch,  Interior  Ballis¬ 
tics  Laboratory,  U.S.  Army  Ballistic  Research  Laboratory,  provided  in-bore 
radar  doppler  data  for  firings  of  two  distinct  projectile  designs  from  a 
specially  configured  37mm  weapon.  Miss  Susan  A.  Coates  of  the  Mechanics  and 
Structures  Branch  accomplished  digital  conversion  of  the  analog  data  at  the 
Laboratory  Experimental  Research  Facility  of  the  Interior  Ballistics  Division 
of  BRL.  Descriptions  of  the  projectile  designs,  test  weapon  and  data  base 
generated  are  presented  in  Section  4.1. 

A  computer  program  was  developed  for  the  specific  purpose  of  model  corre¬ 
lation.  A  description  of  the  program  and  its  required  inputs  is  presented  in 
Section  4.2.  Results  of  the  correlation  effort  are  presented  in  Section  4.3; 
while  several  general  observations  are  presented  in  Section  4.4. 

4.1  Description  of  Test  Weapon,  Ammunition  and  Data  Base 

The  experimental  test  weapon  consisted  of  a  fully  instrumented  37mm  rifled 
gun  tube,  with  a  constant  twist  of  25  calibers  per  turn  and  an  in-bore  projec¬ 
tile  travel  of  72.5  inches.  In  order  to  minimize  gun  tube  motion,  the  system 
design  incorporated  a  heavy,  fixed  collar  supported  near  the  muzzle.  The  test 
weapon  is  shown  in  Figure  7.  The  two  37mm  projectile  designs  consisted  of 
solid,  steel  cores  with  copper  rotating  bands  and  sub-caliber,  steel  bourrelets 
as  '.wn,  along  with  pertinent  data,  in  Figure  8. 

Analog  chamber  pressure  and  in-bore  radar  doppler  data  were  recorded  on 
magnetic  tape  for  all  firings  conducted.  Digital  conversion  and  reduction  were 
accomplished  as  noted  above.  Representative  rounds  for  which  all  required  data 
were  available  were  selected  for  the  purposes  of  model  correlation.  Tabular 
data  for  the  rounds  selected,  consisting  of  pressure,  displacement  and  veloc¬ 
ity-time  histories,  were  transmitted  to  S&D  Dynamics,  Inc. 

The  selected  data  consisted  of  one  record  for  each  projectile  design  whose 
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frequency- time  ("waterfall")  plot  indicated  little  or  no  evidence  of  balloting, 
as  well  as  one  record  for  a  short  ^ound  whose  "waterfall"  plot  indicated  bal¬ 
loting.  As  shown  in  Figure  9,  which  represents  a  portion  of  the  latter  "water¬ 
fall"  plot,  balloting  is  indicated  by  the  presence  of  secondary  peaks  just 
behind  the  main  frequency-ridge  (which  is  directly  proportional  to  the  instan¬ 
taneous  projectile  velocity  along  the  gun  tube  axis).  It  is  noted  that  al¬ 
though  the  secondary  peaks  have  not  as  yet  been  quantified  in  terms  of  projec¬ 
tile  pitch  and  yaw  motions,  it  is  reasonable  to  use  the  time  periods  between 
peaks  as  indications  of  the  frequency  of  balloting.  (Refs.  9,  10) 

4.2  Computer  Program 

For  the  specific  purpose  of  this  correlation,  the  formulation  presented  in 
Section  2  was  programmed  under  the  assumption  of  a  rigid,  straight  gun  tube; 
isolating  projectile  motion  from  gun  tube  motion,  in  accordance  with  the  objec¬ 
tive  of  the  experimental  arrangement.  The  computer  program,  written  in  MBASIC, 
was  debugged  and  executed  on  an  Apple  11+”  computer,  with  CP/M™  operating  sys¬ 
tem.  With  reference  to  Section  2,  the  general  equations  of  motion,  namely 
equations  (35)  and  (36),  were  restructured  into  a  set  of  simultaneous,  non¬ 
linear,  first  order,  ordinary  differential  equations  in  the  projectile  motion 
parameters,  xp,  v  ,  ycg,  yfg,  zcg,  zcg,  tj/,  $,  &,  $,  and  q>Q.  This  set,  with  the 
applied  loads  and  moments  defined  in  Section  2.3,  and  subject  to  the  above 
assumption,  was  integrated  via  a  fourth-order,  fixed  time-step,  Runge-Kutta 
integration  scheme.  Additional  programming  was  introduced  to  detect  the  onset  of 
bourrelet/bore  impact,  to  apply  the  "jump"  conditions  resulting  from  impact, 
and  to  restart  the  integration  scheme  following  each  impact. 

Required  inputs  to  the  program  included  geometric  and  inertia  characteris¬ 
tics  of  the  projectile  design  considered,  as  indicated  in  Figure  8,  as  well  as 
corresponding  tabular  chamber  pressure-time  data.  Elements  of  the  inertia 
matrix  for  each  design  were  hand-calculated  based  on  data  supplied.  In  addi¬ 
tion,  it  was  necessary  to  prescribe  the  rotating-band  "spring"  stiffness  per 
unit  surface  area,  k;  the  coefficients  of  Coulomb  friction  at  the  rotating- 
band/bore  and  bourrelet/bore  interfaces,  namely  y  and  yt  respectively;  the 
instantaneous  projectile  base  pressure,  pb;  and,  finally,  the  "ram"  air  pres¬ 
sure,  p  .  In  the  absence  of  supplied  data  for  these  parameters,  additional 
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modeling  was  undertaken  to  quantify  k,  p.  and  p  ;  whereas,  y  and  y  were  deter- 
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mined  by  parametric  matching  studies. 

In  accordance  with  the  assumptions  introduced  in  Section  2.3.2,  an  elas¬ 
tic/plastic  analysis  of  the  rotating  band,  treating  the  band  as  a  thick-walled 
tube  with  plastic  boundary  initially  at  the  bore  interface  (Ref.  11),  yielded  the 
expression 

[(1  +  v)  r2  +  (1  -  \>)  r?3 

k  =  E  - 2 - —  (83) 

(1  -v2)  r  (r2  -  r.2) 

o  o  1 


where  E  denotes  Young's  modulus  for  the  rotating-band  material,  v  denotes 
Poisson's  ratio,  and  r.  and  r  denote  respectively  the  inner  and  outer  band 
radii.  Furthermore,  it  should  be  noted  that  the  analysis  indicated  only  a 
negligible  increase  in  stiffness  due  to  radially  inward  propagation  of  the 
plastic  zone  for  the  range  of  radial  displacements  anticipated.  Hence,  the 
expression  for  k  as  given  in  equation  (83)  was  taken  as  a  constant  throughout 
the  in-bore  travel. 

The  instantaneous  base  pressure,  p^,  was  prescribed  in  conventional  inte¬ 
rior  ballistics  terms  (Ref.  12,  for  example)  as  a  function  of  the  known  chamber 
pressure,  pc,  in  the  form 


(84) 


where  c  denotes  the  propellant  charge. 

The  "ram"  air  pressure,  p  ,  was  related  to  the  instantaneous  projectile 
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velocity,  Vp,  via  quasi-steady  application  of  the  steady  state,  one-dimensional 
shock  tube  relations  (Ref.  13),  from  which  there  results 


— -  [ (y  +  1 )  M  +  /(y  +  1 )2  M2  +  16  ]  p 
4  P  p  o 


where  pQ  denotes  ambienc  air  pressure,  y  denotes  the  specific  heat  ratio  of  air 
and  Mp  denotes  the  ratio  of  projectile  speed  to  ambient  air  sound  speed. 

Additional  calculated  inputs  and  functional  relations  required  by  the 
program  are  given  in  Figure  10.  Input  to  the  program  is  completed  by  specifi¬ 
cation  of  the  friction  coefficients,  y  and  Pj,  as  defined  above,  and  the  orien¬ 
tation  of  the  projectile  immediately  subsequent  to  full  engraving  (defining  its 
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initial  state) .  It  is  noted  that  non-zero  values  of  projectile  linear  and/or 
angular  displacements  in  its  initial  state  express  the  amount  of  misalignment 
generated  during  the  engraving  process.  Furthermore,  since  this  "malengraving" 
entails  irreversible  plastic  flow,  elastic  restoring  forces  proportional  to 
these  initial  displacements  do  not  exist.  Consequently,  the  expressions  for 
the  rotating-band/bore  interfacial  loads  and  moments  are  modified  within  the 
program  by  subtracting  that  portion  attributable  to  the  initial  state  of  the 
projectile. 

Program  output  included  instantaneous  maximum  radial  displacements  and 
angular  orientations  at  the  rotating  band  and  bourrelet,  as  well  as  linear  and 
angular  displacements  and  velocities  for  all  six  degrees-of-freedom  of  the  pro¬ 
jectile.  In  addition,  bourrelet  motion  relative  to  the  bore  was  monitored, 
with  the  time  of  each  impact  and  the  magnitude  of  each  impulse  output. 

4.3  Results  of  Correlation 

For  the  purpose  of  quantifying  the  value  of  the  rotating-band/bore  inter¬ 
facial  friction  coefficient,  y,  a  series  of  program  runs  were  performed  for  the 
short  round  whose  experimental  record  indicated  little  or  no  evidence  of  ballot¬ 
ing.  Only  the  value  of  y  was  varied  within  these  runs,  all  other  projectile 
parameters  and  initial  conditions  were  held  fixed.  In  addition,  since  the 
experimental  record  did  not  indicate  balloting,  the  projectile  was  assumed  to 
have  negligible  mass  eccentricity  and  misalignment  following  engraving  (i.e., 
initial  linear  and  angular  projectile  displacements  relative  to  the  bore  were 
assumed  to  be  zero).  Output  from  these  runs  which  prescribe  model  predictions 
of  projectile  displacement  and  velocity-time  histories  are  presented,  along 
with  corresponding  experimental  data,  in  Figures  11  and  12,  respectively.  As 
may  readily  be  seen  from  these  figures,  excellent  correlation  was  achieved  for 
a  friction  level  of  y  =  .05.  Hence,  this  value  was  used  for  all  subsequent 
correlations. 

Referring  to  Figures  11  and  12,  it  is  noted  that  model  output  begins 
immediately  subsequent  to  full  engraving.  For  each  round  considered,  required 
initial  values  of  projectile  displacement  and  velocity  were  determined  via 
comparison  of  experimental  data  with  the  known  projectile  geometry  and  free- 
run.  Alignment  of  the  chamber  pressure-time  history  and  determination  of  its 
initial  value  for  each  firing  were  accomplished  by  graphically  matching  the 
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Figure  12  -  Velocity-Time  Correlation  (Short  Round  w/o  Balloting) 
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pressure-time  curve  to  the  projectile  acceleration-time  curve  derived  from  the 
experimental  velocity  record.  This  procedure  was  repeated  for  each  round  con¬ 
sidered. 

The  program  was  exercised  next  for  comparison  with  the  experimental  record 
corresponding  to  the  long  projectile  which  also  evidenced  little  or  no  ballot¬ 
ing.  Model  correlations  with  experimental  data  for  projectile  displacement  and 
velocity-time  histories  are  presented  in  Figures  13  and  14,  respectively.  As 
is  seen,  excellent  theoretical/experimental  agreement  is  obtained  for  this 
round  using  the  friction  level  previously  determined. 

Considering  transverse  and/or  balloting  motion  in  both  of  the  above  cases, 

model  output  indicated  that  (with  negligible  mass  eccentricity  and  malengrav- 
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ing)  the  maximum  radial  displacement  of  the  projectile  is  of  the  order  of  10 
inches  and  the  maximum  yaw  and  pitch  are  of  the  order  of  10-5  milliradians. 
Since  these  values  are  small  compared  to  the  maximum  allowable  values  deter¬ 
mined  by  bourrelet/bore  clearance  and  rotating  band  thickness,  namely,  of  the 
order  of  10-3  inches  and  1  milliradian  respectively,  it  is  reasonable  to  con¬ 
clude  that  the  model  correlates  well  with  experimental  data  regarding  balloting 
motions  for  both  cases  considered. 

Finally,  the  model  was  exercised  for  the  purpose  of  correlation  with  the 
experimental  record  corresponding  to  the  short  round  with  appreciable  evidence 
of  balloting.  The  model  configuration  for  this  correlation  included  a  mass 
eccentricity  of  .02  inches  and  an  initial  malengraved  state  characterized  by 
an  initial  pitch  of  3  milliradians.  However,  it  is  noted  that  (to  within 
practical  limitations)  there  exist  numerous  combinations  of  eccentricity  and/or 
malengraving  which  give  rise  to  balloting.  The  configuration  selected  resulted 
in  appreciable  balloting  with  multiple  bourrelet/bore  impacts  throughout  the 
in-bore  travel,  and  thereby  illustrated  the  full  capabilities  of  the  model. 

For  this  case,  the  coefficient  of  friction  at  the  bourrelet/bore  interface,  Pj , 
was  taken  equal  to  .3  and  the  coefficient  of  restitution  for  impact,  e,  was 
taken  equal  to  1  (rendering  perfectly  elastic  impacts). 

As  for  the  two  previous  rounds  considered,  theoretical  correlation  with 
experimental  projectile  displacement  and  velocity-time  histories  are  shown  for 
this  balloting  round  in  Figures  15  and  16,  respectively.  Once  again  it  is  seen 
that  excellent  theoretical/experimental  agreement  is  achieved  in  terms  of  lon¬ 
gitudinal  projectile  motion. 
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Balloting  motion  as  predicted  by  the  model  for  this  projectile  is  illus¬ 
trated  in  Figure  17;  which  shows  the  time-dependent  variation  of  the  maximum 
radial  displacement  at  the  bourrelet.  In  addition,  the  positions  of  secondary 
peaks  obtained  from  the  experimental  radar  doppler  "waterfall”  plot  for  this 
round  are  also  indicated  in  this  figure.  Equating  these  peaks  with  maximum 
balloting  (as  indicated  by  bourrelet/bore  impacts)  yields  good  agreement  with 
the  dominant  response  frequency  predicted  by  the  model,  especially  in  the  mid¬ 
section  of  the  in-bore  motion.  Variation  of  the  experimentally  derived  fre¬ 
quency  in  the  neighborhood  of  shot-exit,  and  the  poorer  correlation  in  this 
region,  may  reflect  the  effects  of  bore  wear  near  the  muzzle,  as  well  as  dynam¬ 
ic  erosion  of  the  rotating  band  (neither  of  which  are  accounted  for  within  the 
model).  Model  output  further  indicated  that  the  magnitude  and  frequency  of 
balloting  is  dominated  by  projectile  yaw  and  pitch.  This  is  illustrated  in 
Figure  18,  wherein  yaw  versus  pitch  is  presented  relative  to  the  inertial 
reference  frame.  S'. 

4.4  General  Observations 

In  view  of  the  excellent  overall  correlation  with  the  experimental  data, 
several  general  observations  regarding  model  predictions  for  balloting  rounds 
are  not  only  pertinent,  but  have  important  ramifications  regarding  gun  tube 
motion  and  sl'Ot  accuracy.  In  particular,  model  output  for  balloting  rounds 
reveals  high  frequency  content  and  large  magnitudes  in  projectile/bore  inter¬ 
facial  loading,  as  well  as  transverse  linear  and  angular  velocities. 

Referring  to  Figure  17,  three  distinct  frequencies  are  discernible  in  the 
balloting  motion  of  t!  e  projectile  as  predicted  by  the  model.  The  dominant 
frequency  of  approximately  5,000  Hz  is  attributable  to  pitch  and  yaw  motion. 
Superimposed  on  this  dominant  motion  is  a  higher  frequency  of  approximately 
50,000  Hz,  discernible  primarily  at  the  extremes  of  pitching  motion  and 
attributable  to  radial  displacement  of  the  projectile  c.g.  Finally,  a  much 
lower  frequency  of  approximately  70C  Hz  is  discernible  as  amplitude  modulation 
of  the  dominant  mode  and  is  attributable  to  interior  ballistics  pressure  and 
projectile/bore  interfaci.nl  friction  loads. 

Further  insight  into  the  origin  of  the  higher  frequencies  of  balloting 
motion  may  be  gained  by  examination  of  the  expressions  for  the  transverse 
loads  and  moments  generated  at  the  rotating  band/bore  interface  as  given 
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Figure  17  -  Balloting  Correlation  (BRL  Ident.  5) 


respectively  by  equations  (58)  and  (61).  Referring  to  equation  (58),  an  effec¬ 
tive  radial-load  stiffness,  kr,  may  be  defined  as 

k  =  tt  r.  k'  .  (86) 

r  b 

Similarity,  referring  to  equation  (61),  an  effective  ’’foundation"  moment  stiff¬ 
ness,  CM,  may  be  defined  as 
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Based  on  these  stiffnesses,  theoretical  natural  frequencies  for  radial 

(linear)  motion,  defined  as  k^m^,  and  pitch  and/or  yaw  (angular)  motion 

defined  as  CM/I  and/or  C.,/1  ,  were  calculated  for  the  short  round.  These 
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theoretical  frequencies,  namely  45,000  Hz  (linear)  and  4,500  Hz  (angular), 
correlate  well  with  the  frequencies  observed  in  Figure  17.  (Actual  model  out¬ 
put  frequencies  are  expected  to  be  somewhat  higher  since  the  projectile  is  not 
in  a  free-vibration  mode,  its  response  being  modified  by  bore  impacts.) 

It  should  be  noted  that  model  output  for  the  balloting  projectile  indicat¬ 
ed  that  the  rotating-band/bore  interfacial  loads  and  moments  (associated  with 
the  stiffnesses  given  in  equations  (86)  and  (87))  have  maximum  values  of  the 
order  of  10,000  lb  and  1,000  in-lb,  respectively.  In  contrast,  output  for  non¬ 
balloting  projectiles  indicated  interfacial  loads  three  orders  of  magnitude 
lower.  Similarily,  balloting  projectiles  evidence  much  higher  transverse  lin- 
earand  angular  velocities;  approaching  orders  of  10  in/sec  and  10  rad/sec,  re¬ 
spectively.  Such  relatively  high  transverse  velocities  would  be  expected  to 
adversely  affect  projectile  "jump"  at  shot-exit,  and  hence,  overall  shot 
accuracy. 

In  summary,  model  output  indicates  that  balloting  projectiles  experience 
high  frequency,  large  magnitude  interfacial  loads  determined  primarily  by  the 
physical  characteristics  of  the  rotating  band.  These  loads,  which  cannot  be 
accounted  for  in  point-mass  projectile  approximations,  have  the  potential  to 
affect  gun  tube  motion  via  "feedback." 

Since  the  formulation  presented  in  Section  2  also  addressed  application  to 
smoothbore  gun  tubes,  two  additional  program  runs  were  performed  beyond  the 
correlation  study  in  an  attempt  to  contrast  rifled  and  smoothbore  firings.  One 


run  was  performed  for  a  rifled  gun  tube  without  twist;  a  second  run  was  per¬ 
formed  for  a  smoothbore  gun  tube  by  modifying  the  computer  program  to  reflect 
the  additional  degree  of  freedom,  <po>  in  the  equations  of  motion,  as  well  as 
the  required  changes  in  the  impact  equations  (as  noted  in  Section  2).  For  the 
purpose  of  comparison,  both  runs  utilized  the  same  inputs  and  initial  condi¬ 
tions  as  were  used  for  the  short  balloting  round. 

With  the  exception  of  projectile  spin,  comparison  of  the  output  from  these 
two  runs  showed  identical  results  for  the  motion  of  the  projectile.  The  smooth- 
bore  case  exhibited  a  monotonically  increasing  function  for  (pQ,  with  a  maximum 
value  of  5  mrad/sec  at  shot-exit;  whereas,  the  straight  rifling,  of  course, 
constrained  the  projectile  to  zero  spin.  Comparison  of  either  of  these  runs 
with  the  run  with  rifling  twist  showed  that  balloting  induced  by  the  same  ini¬ 
tial  state  of  malengraving  is  markedly  attenuated  by  gyroscopic  stabilization. 
This  is  seen  in  Figure  19,  which  presents  a  comparison  of  the  instantaneous 
pitch  angle  for  spinning  and  non-spinning  (or  smoothbore)  projectiles.  As  is 
seen  in  Figure  19,  there  is  a  marked  increase  in  pitch  magnitude  and  velocity, 
as  well  as  in  the  number  of  bourrelet/bore  impacts,  for  the  non-spinning  projec¬ 
tile.  Output  further  indicated  that  the  magnitude  of  the  impulse  transmitted 
at  bourrelet/bore  impact  was  also  significantly  greater  for  the  non-spinning 
projectile. 

Finally,  a  comparison  of  the  affect  of  balloting  on  the  longitudinal 
motion  of  a  projectile  is  illustrated  in  Figure  20,  wherein  model  output  are 
presented  for  a  perfectly  aligned,  spinning  projectile  (without  bore  impacts 
or  balloting  motion),  as  well  as  spinning  and  non-spinning,  balloting  projec¬ 
tiles.  As  is  seen  from  this  figure,  balloting  adversely  affects  muzzle  veloc¬ 
ity;  propellant  energy  which  would  otherwise  contribute  to  projectile  muzzle 
velocity  in  a  non-balloting  situation  is  expended  in  sustaining  balloting 
motion. 


Figure  19  -  Pitching  Motion  Comparison 


■  ■■■  Spin  Stabilized 
w/o  Balloting 

Spin  Stabilized 
with  Balloting 

— Balloting  w/o  Spin 


J _ i _ I _ _ _ L- 

1  .0  2.0  3.0 

Time  (ms) 


Figure  20  -  Velocity  Comparison 


SECTION  5 

INCORPORATION  OF  MODEL  WITHIN  DYNACODE-G 

Incorporation  of  the  model  herein  developed  within  DYNACODE-G  requires 
consideration  of  the  nature  of  the  equations  of  motion  and  interfacial  loads 
and  moments  presented  in  Section  2,  as  well  as  the  general  observations  noted 
in  Section  4. 

With  respect  to  the  formulation  presented  in  Section  2,  it  is  noted  that 
projectile/gun- tube  motions  are  fundamentally  coupled  through  kinematic,  as 
well  as  loading  constraints.  Kinematic  coupling  is  manifested  by  terms  enter¬ 
ing  the  projectile  equations  of  motion  which  relate  the  instantaneous  orienta¬ 
tion  and  motion  of  the  reference  frame  S  to  the  reference  frame  S'.  Specif i- 

o 

cally,  the  Euler  angles  and  the  coordinates  are  identi¬ 

fied  with  local  gun-tube  slopes,  torque  and  displacements  defined  in  DYNACODE-G. 
Referring  to  the  projectile,  these  quantities  (and  their  time  derivatives)  may 
be  viewed  as  "driving"  forces  entering  the  projectile  equations  of  motion. 
Physically,  these  "driving"  forces  are  transmitted  to  the  projectile  via  load¬ 
ings  acting  at  the  rotating-band/bore  and  bourrelet/bore  interfaces.  Alter¬ 
nately,  referring  to  the  gun  tube,  projectile  "feedback"  to  gun  tube  motion  is 
manifested  through  equal  and  opposite  reaction  loads  and  moments  which  must  be 
applied  to  the  gun  tube  within  DYNACODE-G;  thereby  "driving"  the  gun  tube. 

Hence,  introduction  of  the  model  within  DYNACODE-G  offers  the  ability  to  assess 
the  mutual  effects  of  projectile/gun- tube  interaction  previously  unattainable 
via  a  point-mass  projectile  description. 

Two  additional  points  are  noteworthy.  As  was  seen  in  Section  4,  a  ballot¬ 
ing  projectile  in  an  otherwise  perfectly  symmetric  and  "ideal"  gun  system  has 
the  potential  to  induce  gun  tube  motion.  Conversely,  gun  tube  motion  has  the 
potential  to  induce  balloting  in  an  otherwise  perfectly  aligned  and  symmetric 
projectile. 

To  fully  account  for  the  above  noted  coupling,  incorporation  of  the  model 
within  DYNACODE-G  requires  simultaneous  integration  of  the  differential  equa¬ 
tions  presented  in  Section  2  with  the  equation  set  characterizing  the  general¬ 
ized,  time-dependent  amplitudes  of  the  gun  system  normal  modes  defined  in 
DYNACODE-G.  In  addition,  the  present  sub-routine  which  computes  the  instanta¬ 
neous  loading  on  the  gun  tube  due  to  motion  of  a  point-mass  projectile  is  to 
be  replaced  by  a  routine  incorporating  the  loads  and  moments  developed  in 


As  was  discussed  in  Section  A,  projectile  balloting  is  dominated  by  fre¬ 
quencies  approaching  the  natural  frequencies  of  the  rotating  band  in  angular 
and  radial  modes  respectively.  Because  of  the  relatively  high  values  of  these 
frequencies,  integration  of  the  projectile  equations  of  motion  for  the  parame¬ 
ters  iJj,  $,  y  ,  and  z  required  an  integration  step-size  of  2.5  psec  to 
eg  eg 

achieve  a  stable,  convergent  solution.  In  contrast,  previous  applications  of 
DYNACODE-G  have  demonstrated  stable,  convergent  solutions  with  an  integration 
step-size  of  25  psec  (accommodating  gun  system  natural  frequencies  to  1,000  Hz). 
In  the  interest  of  computational  economy,  it  is  conceivable  to  integrate  pro¬ 
jectile  parameters  at  the  smaller  step-size,  while  utilizing  interpolation  for 
gun  system  parameters  integrated  at  the  coarser  step.  However,  it  is  highly 
likely  that  the  presence  of  projectile  balloting  will  induce  excitation  of 
higher  order  gun  system  modes  in  greater  proportion  than  have  appeired  in  pre¬ 
vious  applications.  Should  this  be  the  case,  convergence  and  stability  of  the 
solution  for  gun  system  parameters  would  require  a  smaller  integration  step- 
size.  In  light  of  these  comments,  it  is  recommended  that  initial  runs  of  the 
combined  code  be  performed  with  the  smallest,  constant  time-step  leading  to 
stable,  convergent  solutions  for  both  projectile  and  gun  system  parameters. 

Exercising  appropriate  caution  in  integrating  the  combined  code  will 
greatly  enhance  our  understanding  of  the  mutual  interaction  between  projectile 
and  gun  system.  For  example,  program  output  will  be  valuable  in  quantifying 
projectile  ''jump”  conditions  at  shot-exit  (relative  to  either  the  gun  tube  or 
an  inertial  reference).  In  addition,  the  combined  code  could  be  used  on  a 
deterministic  basis  to  quantify  shot- to-shot  variations  resulting,  for  example, 
from  variations  in  interior  ballistics  pressure,  projectile  asymmetries  or  the 
degree  of  projectile  malengraving. 
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SECTION  6 

CONCLUSIONS  AND  RECOMMENDATIONS 

A  rigorous  formulation  of  the  general,  six  degree-of-freedom  motion  of  a 
projectile  of  finite  geometry  and  inertia  traveling  in  a  flexible  (rifled  or 
smobthbore)  gun  tube  has  been  developed.  The  formulation  accommodates  projec¬ 
tile  spin,  mass  eccentricity,  projectile/bore  interfacial  friction,  elastic/ 
plastic  deformation  of  the  projectile  rotating  band,  and  balloting — including 
bourrelet  impact  and  rebound  with  the  bore.  The  formulation  has  been  developed 
in  a  form  compatible  with  DYNACODE-G,  such  that  incorporation  within  the  lscter 
provides  the  ability  to  assess  the  mutual  effects  of  projectile/gun- tube  inter¬ 
action  during  in-bore  motion.  The  formulation  has  been  compared  with  Ooher 
projectile  in-bore  motion  formulations  and  descriptions  appearing  in  recent 
literature,  and  has  been  found  to  be  the  most  general  developed  to  date.  In 
addition,  the  formulation  has  been  correlated  with  experimental  in-bore  radar 
doppler  data  provided  by  the  U.S.  Army  Ballistic  Research  Laboratory  for  firings 
of  a  specially  designed  37mm  weapon,  and  has  been  found  to  provide  excellent 
agreement  for  rounds  which  exhibit  balloting,  as  well  as  rounds  which  exhibit 
little  or  no  evidence  of  balloting. 

Referring  to  the  model  output  presented  in  Section  A  for  the  particular 
projectile  designs  therein  considered,  several  observations  regarding  balloting 
are  in  order.  First,  the  onset  of  balloting  is  dictated  by  the  state  of  the 
projectile  when  fully  engraved  (i.e.,  the  total  amount  of  malengraving)  and/or 
the  presence  of  mass  eccentricity.  Second,  balloting  motion  subsequent  to 
engraving,  and  in  the  absence  of  gun  tube  motion  feedback,  is  dictated  primarily 
by  the  design  of  the  rotating  band  and,  in  particular,  the  relative  magnitude 
of  the  "foundation"  moment  stiffness.  Third,  several  distinct  frequencies  are 
discernible  in  the  balloting  motion  of  the  projectiles  considered:  a  dominant 
frequency  of  approximately  5,000  Hz,  dictated  primarily  by  the  response  fre¬ 
quency  of  the  rotating  band  in  the  "foundation"  moment  mode;  a  higher,  super¬ 
imposed  frequency  of  approximately  50,000  Hz,  discernible  at  the  extrema  of 
displacement  in  the  bourrelet-plane  and  dictated  primarily  by  the  radial  re¬ 
sponse  frequency  of  the  rotating  band;  finally,  a  much  lower  frequency  of 
approximately  700  Hz,  discernible  as  amplitude  modulation  of  the  dominant 
frequency  and  primarily  attributable  to  the  interior  ballistics  pressure  and 
projectile/bore  interfacial  friction  loads.  Fourth,  gyroscopic  stabilization 
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(due  to  spin)  has  a  marked  affect  on  the  attenuation  of  balloting,  affecting 
the  loads  transmitted  through  the  rotating  band  and  bourrelet,  as  well  as  the 
linear  and  angular  velocities  of  the  projectile  throughout  its  in-bore  travel 
and,  in  particular,  at  shot  exit.  Fifth,  overall  balloting  motion  as  predicted 
by  the  model  is  of  more  consistent  nature  in  terms  of  frequency  content  than  is 
indicated  in  the  experimental  data  provided,  especially  as  the  projectile  ap¬ 
proaches  the  neighborhood  of  the  muzzle. 

A  final  observation,  although  speculative  at  this  point,  is  nevertheless 
noteworthy.  Gun  tube  motion  predictions  based  on  theoretical  simulations 
generally  contain  considerably  lower  frequency  content  than  do  corresponding 
experimental  data  (e.g.,  accelerometer  data).  Perhaps  this  disparity,  which 
to  date  has  been  viewed  by  many  investigators  as  an  experimental  anomaly, 
arises  as  a  consequence  of  mutual  projectile/gun-tube  interaction.  The  finding 
noted  in  Section  4,  namely,  that  relatively  high  projectile/bore  interfacial 
loads  are  found  to  exist  in  a  balloting  situation,  lends  credence  to  this 
speculation. 

In  view  of  the  demonstrated  merits  of  the  model  herein  developed  and  its 
potentially  broad  range  of  applicability,  as  well  as  its  potential  for  resolving 
the  theoretical/experimental  disparity  noted  above,  it  is  strongly  recommended 
that  this  study  be  extended  for  the  purpose  of  incorporating  the  model  within 
DYNACODE-G,  as  discussed  in  Section  5.  In  addition,  in  view  of  the  importance 
of  accurately  defining  the  state  of  the  projectile  immediately  following  engrav¬ 
ing,  as  well  as  defining  variations  in  projectile/bore  clearances  due  to  wear 
and  erosion  of  both  the  bore  and  rotating  band  (particularly  in  the  latter 
stages  of  in-bore  travel),  it  is  recommended  that  additional  effort  be  expended 
in  modeling  these  processes  (via  introduction  of  stochastic  techniques,  where 
applicable) . 
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APPENDIX  A 


INSTANTANEOUS  ANGULAR  ORIENTATION 
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APPENDIX  A 

angular  orientation  and  velocity  of  the  intermediate 

AAA 

,  with  coordinates  (x  ,y  ,z  )  and  unit  triad  (i  ,j  ,k  ), 

o  o  o  o  o  o 

ve  to  the  inertial  reference  frame.  S' ,  with  coordinates 

AAA 

triad  ( i' ,j' ,k' ) ,  in  terms  of  the  Euler  angles 

-1. 


Sq  initially  coincident  with  S' ,  the  instantaneous  anguh 


elative  to  S'  at  any  later  time  is  (referring  to  Figure  A- 


ing  S°  to  the  following  consecutive  rotations: 


t  z' ,  bringing  x'  to  its  final  elevation,  and 
ts  intermediate  orientation,  T); 


c  r|,  bringing  £  to  its  final  azimuth,  xq,  and 


.ts  final  azimuth ,  £; 


rt  xq,  bringing  n  to  its  final  orientation,  yQ, 


o  its  final  elevation ,  z  . 

o 

cosines,  i..,  definino  the  transformation  between  S  and 
ij  o 


re  obtained  by  noting  the  relation  between  the  unit  triad. 
a-1  subsequent  to  each  consecutive  rotation,  as  follows: 


tation  .  there  results 
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Eliminating  the  unit  triad  ( i^,i  ,i  }  associated  with  the  intermediate 

s  n  s 

coordinates  (f,,r),£)  there  results  the  transformation  as  given  in  equation  (5), 
with  the  direction  cosines  as  defined  in  equation  (6). 

Referring  again  to  Figure  A-l,  it  is  seen  that  the  angular  velocity  of  Sq 

relative  to  S' ,  namely  Oi  is  given  as 


W  =  4>  k' 
o  o 


A  • 


A 

Jo  ‘ 


(A-7) 


Once  again  using  the  above  transformations  between  unit  triads,  there 
results  in  the  form  as  given  in  equation  (12);  with  components  as  defined 
in  equation  (13). 


NOMENCLATURE 


projected  areas  of  rotating  band  and  bourrelet 
total  translational  acceleration  of  projectile  c.g. 
relative  to  S' 

acceleration  of  projectile  c.g.  relative  to  SQ 

acceleration  of  origin  of  Sq  relative  to  S' 

projectile  acceleration  along  gun  tube  axis 

components  of  a  relative  to  S 
o  o 

effective  "foundation"  moment  stiffness 
propellant  charge 

resultant  incremental  contact  load  acting  on  a  differential 
element  of  the  rotating  band 
components  of  dFc  relative  to  SQ 

Young's  modulus  for  rotating  band  material 

coefficient  of  restitution 

resultant  applied  load  acting  on  projectile 

impulsive  load  at  bourrelet/bore  interface 
"ram"  air  pressure  load 

projectile  base  pressure  load 

resultant  pressure  load 

projectile  weight  load 

components  of  F  relative  to  Sq 

components  of  rotating-band/bore  interfacial  loads 
relative  to  S 

o 

components  of  Fp  relative  to  SQ 
components  of  Fw  relative  to  SQ 
gravitational  acceleration 

I  PREVIOUS  PAGE  9k 
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H  -  total  angular  momentum  of  projectile 

h  -  rotating  band  width 

I  ,  I  .....  I  -  elements  of  projectile  inertia  tensor  relative  to  S 
xx’  yy’  '  yz  H  J 


AAA 

(i'J'.k') 

AAA 

(1o>Jo>kol 


k' 


M 


^  (i, j=1 ,2,3) 

S,°j  ( i , j  =  1  ,2,3) 

M 

M 


M 


-  unit  triad  associated  with  S 

-  unit  triad  associated  with  3’ 

-  unit  triad  associated  with  S 

o 

-  radial  unit-vector  at  bore  surface 

-  elastic  "spring"  stiffness  per  unit  surface  area  of 
rotating  band 

-  "spring"  stiffness  per  unit  circumferential  length  of 
rotating  band 

-  effective  radial-load  stiffness  of  rotating  band 

-  moment  arm  from  projectile  c.g.  to  "ram"  air  load 

-  perpendicular  distance  from  projectile  c.g.  to  rear  face 
of  rotating  band 

-  moment  arm  from  projectile  c.g.  to  base  pressure  load 

-  moment  arm  of  incremental  contact  load  at  rotating-band/ 
bore  interface 

-  perpendicular  distance  from  projectile  c.g.  to  forward 
face  of  bourrelet 

-  direction  cosines  between  S  and  S 

o 

-  direction  cosines  between  S  and  S' 

o 

-  resultant  applied  moment  acting  about  projectile  c.g. 

-  ratio  of  projectile  speed  to  ambient  air  sound  speed 

-  moment  due  to  pressure  loads 


-  components  of  M  relative  to  S 

-  components  relative  to  S  of  moment  due  to  rotating-band/ 
bore  interfacial  contact 

-  components  of  relative  to  S 

-  components  relative  to  SQ  of  moment  due  to  rotating-band/ 
bore  interfacial  contact 


ICC  <03  It-  £,  it,  t,  t,  it,  it,  CO  CO  CO  B  ih 


k-  r- am-*  r  u  i  •„*»  "j" -  a*v  *j<i  — ^ 


'  >i rjv  ^  wt.  •T' «  tf'.;  **-,£  ir*  i'i-  i.ir^min.1  Bturvijr^  rj  ^’H-irciwasw  ■; 


m_ 


1 


« 


-  projectile  mass 

-  "ram"  air  pressure 

-  projectile  base  pressure 

-  chamber  pressure 

-  ambient  air  pressure 

-  radial  load  per  unit  surface  area  of  rotating  band 

-  radially  directed  impulsive  load  acting  at  bourrelet/bore 
interface 

-  position  vector  of  projectile  c.g.  relative  to  origin  of  S’ 

-  bore  radius 

-  position  vector  of  projectile  c.g.  relative  to  origin  of  SQ 

-  inner  radius  of  rotating  band 

-  outer  radius  of  rotating  band 

-  position  vector  of  origin  of  SQ  relative  to  origin  of  S' 

-  position  vector  from  projectile  c.g.  to  point  of  bourrelet/ 
bore  impact 

-  projectile  body-fixed  reference  frame 

-  inertial  reference  frame 

-  intermediate  (gun  tube)  reference  frame 

-  distance  along  gun  tube  axis 

-  rifling  torque-load  per  unit  surface  area  of  rotating  band 

-  time 

-  rifling  twist 

-  total  translational  velocity  of  projectile  c.g.  relative 
to  S' 

-  velocity  of  bourrelet  impact-point  relative  to  the  bore 

-  velocity  of  projectile  c.g.  relative  to  SQ 

-  velocity  of  projectile  c.g.  relative  to  the  gun  tube 

-  velocity  of  origin  of  SQ  relative  to  S’ 


31 


.  .  *  r-  u  r  .j  .  1  ».  1  «j  '.-W  ~w  ' 


CTiii  '  v  -  v  -*->.-»  .-  s  r  J-J  r.  <rj,  .  *  .  -.  VAT? 


V  ,  V  ,  V 

x  *  y  *  z 
o  Jo  o 


(x,y,z) 

(xf ,y’ ,z* ) 


(W*o> 


ycg’  zcg 


Greek  Symbols 
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-  projectile  velocity  along  gun  tube  axis 


components  of  vQ  relative  to  SQ 


orthogonal  coordinates  of  S 
orthogonal  coordinates  of  S' 
orthogonal  coordinates  of  S 


projectile  c.g.  displacements  in  jQ-ko  plane 


-  angular  orientation  of  friction  load  at  rotating-band/ 


bore  interface 


specific  heat  ratio  of  air 

transverse  displacement  of  any  point  in  a  plane 
perpendicular  to  projectile  geometric  axis 
maximum  radial  displacement 


-  initial  radial  compression  of  rotating  band 


-  radial  component  of  6 


projectile  c.g.  offset  from  geometric  axis 
projectile  yaw  relative  to  S’ 

friction  coefficient  at  rotating-band/bore  interface 
friction  coefficient  at  bourrelet/bore  interface 


-  Poisson's  ratio 


distance  from  projectile  c.g.  to  any  plane  perpendicular 
to  geometric  axis  of  projectile 

impulsive  torque  generated  at  rotating-band/bore  interface 
due  to  bourrelet/bore  impact 


Euler  angle  of  S  relative  to  SQ;  also,  angle  between 


projection  of  yQ-axis  and  line  from  gun  tube  centerline 


to  point  in  plane  of  interest 

angular  orientation  of  maximum  radial  displacement 


Euler  angle  of  SQ  relative  to  S';  also,  projectile 


spin  angle 


-  angular  orientation  of  bourrelet/bore  impact-point 
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-  Euler  angles  of  S  relative  to  SQ 

-  projectile  pitch  relative  to  S* 

-  Euler  angles  of  Sq  relative  to  S' 

-  total  angular  velocity  of  projectile  relative  to  S'  • 

-  total  angular  acceleration  of  projectile  relative  to  S' 

-  angular  velocity  of  projectile  relative  to  SQ 

-  angular  acceleration  of  projectile  relative  to  SQ 

-  angular  velocity  of  SQ  relative  to  S' 
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-  angular  acceleration  of  Sq  relative  to  S' 

-  components  of  w  relative  to  S 

-  components  of  to  relative  to  S 

-  components  of  to  relative  to  S 

eg 

• 

-  components  of  w  relative  to  S 

eg 

-  components  of  uio  relative  to  SQ 


w  ,  to  to  -  components  of  w  relative  to  S 

x  y  z  o  o 
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Additional  Notation 

A 

(  )  -  unit  vector;  also,  impulsive  load 

(  )  -  vector  quantity 

3(  )  -  partial  derivative 

d(  )  -  total  (convective)  derivative 

(  )  -  derivative  with  respect  to  time 

A(  )  -  "jump"  (sudden  increment)  in  parameter  due  to 

bourrelet/bore  impact 

(  )+  -  value  of  parameter  immediately  subsequent  to  impact 

(  )”  -  value  of  parameter  just  prior  to  impact 
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